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than  predicted  by  linear  theory.  This  report  examines  some  possible  explanations  of 
that  observation.  In  particular,  we  study  the  effect  of  non-linearities  in  the 
equation  of  state  due  to  entropy  generation,  the  variation  of  sound  generating 
efficiency  with  time  for  subsonic  and  transonic  sources,  and  the  validity  of  linear 
superposition  for  a  moving  thermoacoustic  source. 

Our  analytical  study  of  the  effect  of  a  nonlinear  equation  of  state  on  sound 
generation  for  a  transonic  source  shows  that  such  considerations  are  negligible, 
being  several  orders  of  magnitude  smaller  than  the  linear  driving  term.  Numerical 
results  for  the  efficiency  of  the  sound  generating  process  as  a  function  of  time 
indicate  that  the  efficiency  is  smaller  for  a  transonic  source  than  for  a  subsonic 
source.  Hence  there  is  no  support  for  a  pumping  principle.  Finally,  experimental 
work  shows  that  the  superposition  principle  holds  better  for  a  transonic  source  than 
for  either  a  subsonic  or  a  supersonic  source.  This  surprising  result  rules  out 
nonlinearities  as  the  cause  of  the  observed  boost  in  sound  pressure  level  for  the 
transonic  source. 

•In  Part  II  of  this  report  we  discuss  how  the  sound  signal  generated  by  a  laser 
source  is  affected  by  various  source  and  receiver  parameters/-"  Previously,  work  on 
the  laser  source  has  been  limited  to  determination  of  the  root  mean  square  sound 
pressure  level  as  a  function  of  the  angle  between  the  receiver  and  the  initial 
source  location.  Since  the  source  is  moved,  the  signal  received  by  the  hydrophone 
varies  markedly  as  the  source  and  geometrical  parameters  of  the  systems  vary.  This 
study  determines  the  relative  importance  of  these  parameters.  The  results  could 
prove  to  be  very  useful  if  the  thermoacoustic  source  is  ever  employed  as  an 
underwater  communications  device. 

In  tnis  study  a  numerical  code  was  used  to  simulate  the  acoustic  signal  observed  by 
a  receiver  as  the  source  parameters  (such  as  velocity  and  modulation  frequency)  and 
geometrical  parameters  (such  as  the  initial  range  and  angle  between  the  source  and 
receiver)  are  varied--  Temporal  and  amplitude  signal  properties  are  examined.  It  is 
shown  that  the  amplitude  properties  are,  in  general,  more  sensitive  to  the  changes 
in  the  system  than  the  temporal  properties.  Further,  it  is  demonstrated  that  the 
source  speed  arid  the  angle  between  the  source  and  hydrophone  are  the  most  critical 
parameters . 
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Part  I 


SOUND  GENERATION  BY  A 
TRANSONIC  ALLY  MOVING 
LASER  SOURCE 


CHAPTER  1 


INTRODUCTION 


Laser  induced  sound  has  been  the  subject  of  research  around  the  world 
for  the  past  20  years.  Much  of  this  work  is  summarized  in  two  review  articles  by 
Lvamshev  and  Sedov1  and  Lyamshev  and  Naugol’nykh.2  There  are  five  mecha¬ 
nisms  for  the  generation  of  sound  by  a  laser:  electrostriction,  thermal  expansion, 
surface  evaporation,  explosive  boiling,  and  optical  breakdown.  Which  of  these 
sound  generating  processes  is  dominant  depends  on  the  laser  and  fluid  parame¬ 
ters,  especially  the  laser  beam  energy  density.  Although  the  optoacoustic  energy 
conversion  efficiency  of  the  mechanism  of  thermal  expansion  is  relatively  low, 
this  is  the  mechanism  that  has  been  primarily  investigated  because  the  process  is 
easily  controllable,  and  is  capable  of  producing  a  highly  directional  sound  source. 

In  the  thermal  mechanism,  the  heating  of  the  medium  by  a  laser  causes 
a  temperature  rise  and  a  corresponding  density  drop.  Modulation  of  the  laser  at 
frequency  /  thus  generates  sound  waves  at  that  frequency. 

The  first  substantial  work  on  the  generation  of  sound  by  heat  may  be 
attributed  to  Ingard3  who  derived  a  wave  equation  describing  the  sound  field  of 
a  lossless  medium  containing  heat  sources.  Application  of  this  work  to  the  sound 
generated  by  a  laser  was  first  described  by  Westervelt  and  Larson,4  who  derived 
the  equations  which  describe  the  process.  Muir,  Culbertson,  and  ('lynch''  later 


conducted  the  first  experiments  to  verify  the  predictions  of  Westervelt  and  Larson. 
This  experimental  work  was  limited  to  the  case  where  the  laser  beam  is  stationary 
so  that  heat  deposition  is  confined  to  a  fixed  portion  of  the  medium.  The  energy 
conversion  efficiency  for  this  scheme  was  found  to  be  quite  low,  prompting  further 
study  into  means  of  enhancing  the  sound  generation  efficiency.  Bunkin  et  al.,6’' 
Bozhkov  et  al.,8  and  Lyamshev  and  Sedov9  suggested  that  a  means  u'  increasing 
the  efficiency  of  sound  generation  is  through  the  use  of  a  moving  thermoacoustic 
source,  in  which  the  laser  beam  is  scanned  across  a  volume  of  the  liquid.  They 
argued  that  moving  the  source  at  Mach  1  could  boost  the  efficiency  of  the  process 
by  causing  all  of  the  sound  pressures  in  the  direction  of  motion  to  add  coherently. 

Berthelot  and  Busch- Vishniac10-12  conducted  an  extensive  analytical 
and  experimental  study  of  moving  thermoacoustic  sources.  Unlike  the  analy¬ 
ses  reported  in  the  Soviet  literature,  their  work  is  valid  for  both  the  nearfield  and 
the  farfield,  includes  the  effect  of  an  air-water  interface,  uses  a  time  domain  ap¬ 
proach,  and  includes  the  effects  of  a  finite  beamwidth.  Their  theoretical  model  for 
moving  thermoacoustic  sources  is  based  on  the  linear  wave  equation  derived  by 
Larson.13  Their  approach  uses  numerical  predictions  based  on  a  convolution-type 
summation  between  the  impulse  response  of  the  system  and  the  thermoacous¬ 
tic  source  strength.  This  theory  is  well  supported  by  experimental  results  for 
the  cases  of  subsonic  and  supersonic  motions  of  the  source.  However,  a  ma¬ 
jor  discrepancy  was  found  for  the  transonic  case  where  the  theoretical  model 
significantly  underestimates  the  measured  sound  pressure  level.  Figure  (1-1), 
reproduced  from  Ref.  12,  is  typical  of  the  results  obtained  by  Berthelot  and 
Busch- Vishniac.  Note  that  there  is  roughly  a  25  dB  disagreement  between  the 
predicted  and  measured  results  near  Mach  1.  Based  on  these  results  and  articles 
published  in  Soviet  journals211-10  it  was  suspected  that  nonlinear  effects  became 
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important  near  Mach  1.  In  particular,  it  was  suggested  that  the  coefficient  of 
nonlinearity  and  the  specific  heat  of  the  medium  may  both  change  as  the  source 
moves,  due  to  the  associated  temperature  changes. 

In  recent  years,  several  theories  have  been  proposed  for  increasing  the 
transduction  efficiencies  of  the  thermoacoustic  generation  of  sound.  One  ot  them, 
discussed  by  Hsieh17  and  Pierce  and  Hsieh,18  predicted  that  the  rate  of  acoustic- 
energy  production  per  unit  volume  is  proportional  to  the  product  of  the  rate  of 
heat  addition  per  unit  volume  and  the  local  acoustic  pressure  presently  existing 
at  the  same  point.  According  to  this  principle,  for  maximum  efficiency  the  laser 
energy  deposition  scheme  should  create  a  traveling  acoustic  wave  and  then  sys¬ 
tematically  pump  up  its  amplitude  by  always  adding  energy  at  places  where  the 
acoustic  pressure  has  a  local  maximum. 

This  report  discusses  an  attempt  to  determine  whether  nonlinear  inter¬ 
actions  may  explain  the  observed  discrepancy  near  Mach  1.  The  wwk  is  partly 
analytical,  partly  numerical,  and  partly  experimental.  The  numerical  work  is 
based  on  the  work  by  Berthelot.19 

Chapter  2  examines  the  effect  of  including  fluid  properties  which  vary 
with  temperature  through  a  nonlinear  equation  of  state.  A  linear,  second  order 
wave  equation  is  derived  which  contains  terms  to  second  order  in  the  specific  heat 
and  coefficient  of  thermal  expansion  of  the  medium.  The  effects  of  these  terms 
in  the  wave  equation  are  examined  theoretically  and  numerically.  Pressure-time 
waveforms,  generated  with  and  without  these  terms,  are  compared. 

The  efficiency  of  thermoacoustic  radiation  is  the  subject  of  Chapter  3. 
The  expression  for  efficiency  defined  by  Pierce20  is  evaluated  numerically  as  a 
function  of  time  for  subsonic  and  transonic  motions  of  the  thermoacoustic  source 
on  the  water  surface. 


Chapter  4  presents  experimental  observations  designed  to  determine  the 
importance  of  nonlinear  effects.  Verification  of  the  superposition  principle  is  dis¬ 
cussed  for  sources  moving  at  Mach  numbers  less  than,  equal  to,  and  greater  than 
one.  A  detailed  description  of  the  experiment  and  equipment  used  is  also  given. 

Finally,  the  results  of  this  part  of  the  study  are  summarized  and  relevant 
conclusions  are  drawn  in  Chapter  5. 


CHAPTER  2 


DERIVATION  AND  ANALYSIS  OF  THE  LINEAR  SECOND 
ORDER  WAVE  EQUATION  FOR  THERMOACOUSTIC  SOUND 

GENERATION 


As  mentioned  earlier,  the  theoretical  model  developed  by  Berthelot  and 
Busch-Vishniac12  for  moving  thermoacoustic  sources  underestimates  the  mea¬ 
sured  sound  pressure  level  when  the  source  moves  at  speeds  near  Mach  1.  Several 
attempts  have  been  made  to  explain  the  discrepancy  between  theoretical  and  ex¬ 
perimental  results  for  transonic  motion  of  the  source  on  the  water  surface.  One 
such  theory,  proposed  by  Soviet  scientists, 214-16  put  forward  the  idea  that  the 
thermodynamic  variables  of  the  medium  cannot  be  taken  as  invariant.  They  sug¬ 
gested  that  the  time  rate  of  change  of  these  properties  becomes  too  important  to 
be  neglected,  especially  for  source  speeds  near  Mach  1.  In  particular,  concern  was 
expressed  regarding  the  specific  heat,  Cp,  and  coefficient  of  thermal  expansion.  3. 

In  this  chapter  we  examine  the  possible  effects  of  a  temperature  (and 
time)  varying  set  of  properties  of  the  medium.  We  develop  the  partial  differential 
equation  describing  sound  generation  in  the  medium  from  the  standard,  linear 
continuity  and  conservation  of  momentum  equations,  and  a  nonlinear  state  equa 
tion.  I  he  state  equation  we  use  includes  the  normal  first  order  terms  plus  the 
second  order  term  which  goes  as  the  square  of  the  entropy  generated  due  to  the 


heat  input  to  the  system.  All  other  second  order  terms  in  the  equation  of  state 
are  neglected.  In  this  manner  we  specifically  examine  only  the  nonlinear  effect 
attributable  to  the  laser  source  supplying  sufficient  heat  to  cause  a  change  in  the 
properties  of  the  medium.  The  order  of  magnitude  of  each  of  the  resulting  second 
order  source  terms  in  the  acoustic  generation  equation  is  estimated  to  find  those 
which  dominate.  I  he  pressure-time  response  is  computed  and  plotted  for  each  of 
the  predominant  source  terms,  and  compared  to  the  first  order  source  term. 

2-1  DERIVATION  OF  THE  ENERGY  AND  WAVE  EQUATIONS 

Let  P  =  (po  +  p)  be  the  total  pressure  in  the  medium,  where  p0  is  the 
static  pressure  and  p  is  the  acoustic  pressure  generated  by  the  thermal  mecha¬ 
nism.  Similarly  let  p  =  (p0  +  8p)  be  the  total  density,  T  =  (T0  - f  ST)  the  total 
temperature,  and  s  =  (s0  +  <$s)  the  total  specific  entropy  of  the  fluid  medium. 
The  terms  p0,  To,  and  So  refer  to  the  equilibrium  properties  of  the  medium  and 
Sp.  ST,  and  Ss  are  those  portions  resulting  from  acoustic  disturbances. 

The  linearized  continuity  and  momentum  equations  are 

Pt  T  Po ^  •  u  -*  0  (2-1.1) 

PoUt  +  Vp  =  0  (2-1.2) 

Here  u  is  the  particle  velocity  in  the  medium  and  the  subscript  t  denotes  a 
time  derivative.  For  simplicity,  it  has  been  assumed  that  the  medium  is  lossless, 
nonconvecting.  uniform,  and  irrotational. 

In  the  standard  equation  of  state  used  in  acoustics,  one  assumes  that 
the  density  of  the  medium  may  be  expressed  solely  in  terms  of  the  pressure.  For 
thermally  generated  sound  it  is  necessary  to  recognize  that  the  density  of  the 
medium  is.  in  general,  a  function  of  two  independent  thermodynamic  variables 


9 


S' 
■  ^ 


rather  than  one.  We  chose  to  allow  the  independent  variables  to  be  the  specific 
entropy  and  pressure.  Hence  the  state  equation  for  this  thermal  process  is 


P  =  P(P-*)  ■ 

We  use  a  second  order  Taylor  series  expansion  to  approximate  6p  by 


r  VP 

SP=dP 


dp  I  1  d2p 

p+d;\Js  +  2d? 


(6s 


(2-1.3) 


(2-1.4) 


where  we  have  neglected  the  other  second  order  terms  because  we  assume  that 
variations  in  s  resulting  from  heating  of  the  medium  dominate  over  variations 
in  p.  By  definition,  the  sound  speed,  c0,  is  given  by  the  relation 
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dp 


Therefore. 


djp_ 

dP\ 


0  co 


he  coefficient  of  thermal  expansion,  /?,  is  defined  as 


p  9T 


or 


,?=  _i  dp 

p  ds 


ds 

df 


But,  from  thermodynamics, 


c. 


ds 


-JL  = 

r  dr  I, 


where  Cp  is  the  specific  heat  of  the  medium.  Therefore, 

?  dp\  £p_ 

p  ds  \p  T 


and 
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(2-1.5) 


(2-1.6) 


(2-1.7) 


(2-1.8) 


(2-1.9) 
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Differentiating  Eq.  (2-  1.1 1)  with  respect  to  s ,  we  get 
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(2  1-12) 


We  may  use  Eq.  (2  1.11)  to  express  in  Eq.  (2-1.12),  but  P-l  .  % 


aa  I 
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ds  i 


and  '-Tf-  are,  as  yet,  unknown.  In  general,  T,  /?,  and  Cp  are  functions  of  both 
P  and  s.  However,  since  we  are  assuming  that  the  entropy  production  dominates 
over  sound  generation,  and  since  these  coefficients  appear  in  the  highest  retained 
order  of  our  approximation  for  6p,  we  may  approximate  T.  /3,  and  Cp  as  functions 
of  s  only.  Therefore  we  may  write  |^|  6s  as  63  and  similarly  define  6T  and  6CP. 
Then,  our  expansion  for  6p  produces 
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(2-1.13) 


where  §^Jo  6s  is  approximated  to  first  order  as  6p- Differentiating  Eq.  (2-1.13) 
with  respect  to  time  yields 
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(2-1-14) 


Eq.  (2-1.14)  contains  terms  which  explicitly  depend  upon  the  time  rate 
of  change  of  entropy  produced  by  the  laser  heating  of  the  medium.  We  prefer 
to  express  the  density  perturbation  instead  in  terms  of  the  heat  supplied  by  the 
laser,  so  we  use  the  second  law  of  thermodynamics  to  relate  the  entropy  and  heat. 
Note  that  this  assumes  that  all  of  the  heat  results  in  an  entropy  change. 

q  a  1 

at  =  ~r  =  l - 71  T7r~T TFT  •  (2-1.15) 

pd  (Po  +  6p)  (T0  +  61 ) 

By  taking  the  binomial  expansion  of  the  denominator  and  retaining  the  leading 
terms,  we  have 
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(2  1.16) 


where  q  is  the  heat  input  per  unit  volume  per  unit  time.  Substituting  Eq.  (2  1.16) 
into  Eq.  (2-1.14)  and  keeping  leading  terms,  Eq.  (2  1.14)  simplifies  to 
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which  may  be  rewritten  as 
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r o  ^ou2)  +  ^  •  (pu) = ~—pt  . 

V  2  / 1  po 


(2  1.17) 


Multiplying  Eq.  (2-1.1)  throughout  by  p,  we  have 

P(V  •  u)  =  -  —pt  • 

Po 

Taking  u-  of  Eq.  (2-1.2)  yields 

u  •  (p0u()  =  -u  ■  Vp  =  -V  •  (pu)  +  p(V  •  u)  . 

Therefore 

p(V  •  u)  =  u  •  (p0u,)  +  V  •  (pu)  . 
Equating  p(V  -u)  from  Eqs.  (2-1.18)  and  (2-1.20), 
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(2-1.18) 


(2-1.19) 


(2-1.20) 


(2-1.21) 


(2-  1.22) 


Equation  (2-1.22)  may  be  interpreted  as  an  energy  equation.  The  terms  on  the 
left  hand  side  are  the  kinetic  and  potential  energies  per  unit  volume  per  unit 
time,  and  the  term  on  the  right  hand  side  is  the  source  factor.  Note  that  since 
the  source  term  depends  upon  rlio,.  the  second  order  terms  in  the  equation  of 
state  have  the  affect  of  acting  as  new  source  terms. 
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Substituting  for  p,  from  Eq.  (2-1.17)  in  Eq.  (2-1.22)  yields 
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(2  1.23) 


It  is  this  equation,  combined  with  the  conservation  of  linear  momentum  and 
conservation  of  mass,  that  will  be  used  to  produce  a  wave  equation  for  the  thermo¬ 
acoustic  source.  This  energy  equation  is  also  useful  in  discussing  the  efficiency  of 
thermoacoustic  generation  as  will  be  discussed  in  the  next  chapter.  For  now,  it  is 
sufficient  to  note  that  the  leading  order  source  term  in  Eq.  (2-1.23)  goes  as  the 
product  of  q  and  p.  and  that  all  other  source  terms  involve  the  product  of  p  and 
other  thermodynamic  variables. 

Consider  the  energy  equation,  Eq.  (2-1.23).  Multiplying  throughout  by 
(—po/p)  and  regrouping  produces 
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Using  the  momentum  equation,  Eq.  (2-1.2),  Eq.  (2-1.24)  reduces  to 
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Differentiating  this  equation  with  respect  to  f  yields 
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Equation  (2  -1.26)  may  be  reduced  using  the  divergence  of  the  momen¬ 
tum  equation  (Eq.  (2-1.2))  to  yield  the  wave  equation  including  second  order 
terms  attributed  to  time  varying  thermodynamic  parameters: 
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Note  that  although  the  left  hand  side  of  this  equation  contains  the  stan¬ 
dard  acoustic  terms,  there  are  also  pressure  dependent  terms  on  the  right  hand 
side  of  Eq.  (2-1.27).  Note  also  that  there  are  20  new  terms  which  appear  on 
the  source  side  of  the  equation  and  which  result  directly  from  the  inclusion  of  a 
nonlinear  equation  of  state.  The  remaining  term  on  the  right  hand  side  is  the 
leading  order  source  term  which  appears  in  the  West erveld- Larson  equation. 

The  wave  equation  derived  above  is  too  complicated  to  be  analyzed 
theoretically  or  numerically.  Its  complexity  can  be  reduced  by  making  certain 
simplifying  assumptions,  and  examining  the  orders  of  magnitude  of  some  of  the 
source  terms.  There  are  two  operations  which  may  be  readily  applied  to  reduce 
t he  complexity  of  Eq.  (2  1 .27):  eliminat ing  p  using  t he  equat ion  of  state,  and  ex¬ 
pressing  all  time  derivatives  of  thermodynamic  quantities  in  terms  of  derivatives 
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with  respect  to  temperature  and  time  derivatives  of  temperature. 

To  eliminate  p  from  Eq.  (2  1.27),  we  rewrite  dp.  pt.  and  pn  using  the 
leading  order  terms  in  the  state  equation.  (This  is  permissible  since  the  terms 
involving  rho  are  all  second  order.)  Making  this  substitution  in  Eq.  (2  1.27) 
vields 
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(2-1.28) 


where  we  recognize  q  =  p7’.s,  &  poTo$t- 

If  the  thermodynamic  quantities  Cp  and  3  are  functions  of  a  single  ther¬ 
modynamic  variable,  as  assumed,  then  we  may  express  them  as  CP(T)  and  3(1'). 
Using  the  chain  rule,  then,  we  may  rewrite  Eq.  (2  -1.28)  as 
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where  we  have  retained  only  terms  to  second  order.  It  is  this  simplified  wave 
equation  whirl)  we  will  examine  in  order  to  determine  the  relative  magnitudes  of 
the  various  source  terms. 


2-2  MAGNITUDES  OF  THE  SOURCE  TERMS  IN  THE  WAVE 
EQUATION 

In  order  to  compare  the  magnitudes  of  the  source  terms  appearing  in 
Eq.  (2-1.29)  it  is  convenient  to  relate  the  perturbations  in  the  thermodynamic 
variables  to  q.  The  expression  for  q  in  terms  of  the  intensity  (I)  of  the  laser  beam 
has  been  shown  bv  Larson13  to  be 


q=  -V  ■  ( I{x,y,z,t )) 


(2-2.1) 


Considering  the  exponential  variation  of  intensity  with  depth,  to  be  the  only 
spatial  variation  of  intensity  within  the  laser  beam,  we  have 


q  =  [/t/(f)exp(-az)]  S8(x  -  xc)6{y  -  ye)  , 
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(2-2.2) 


where  A  is  a  constant  for  a  given  laser,  S  is  the  cross-sectional  area  of  the  beam, 
q  is  the  optical  absorption  coefficient,  and  8(x  —  xc)  and  6(y  —  yc)  are  Dirac  delta 
functions.  From  the  above  equation,  for  a  given  a  and  at  a  given  depth, 


q  oc  I(t)  . 


For  the  laser  used  in  our  research 


I{t)  oc  —  exp(— 5(</ rp))  , 
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where  rp  is  the  laser  pulse  duration.  Therefore 


q  oc  —  exp(-5(f/rp))  . 


(2-2.3) 
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I  he  above  expressions  allow  for  q  and  r/,  to  he  expressed  in  terms  of 
known  quantities.  However,  we  may  also  write  tj  as 
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Equation  (2  2.6)  may  be  approximated  to  leading  order  by 
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which  may  be  solved  for  Tt. 
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Integrating  this  equation  with  respect  to  time  yields  an  expression  for  the  temper¬ 
ature  rise  of  a  volume  caused  by  the  laser  heating. 
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Hence,  we  may  express  T,  and  ST  in  terms  of  known  quantities. 

Since  Cp  and  i3  are  assumed  to  be  only  functions  of  T,  6CP  and  SO  can 
be  expressed  as 
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and 
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With  these  definitions  expressions  have  berm  obtained  for  all  of  the  source  terms 
in  the  simplified  wave  equation,  except  for  those  involving  the  acoustic  pressure. 
Since  Eq.  (2  1.29)  is  linear  in  pressure,  we  estimate  the  strength  of  the  pressure 
dependent  source  terms  by  first  determining  the  pressure  response  to  the  leading 
order  source  term,  and  then  using  this  pressure  in  the  second  order  source  terms. 
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Although  this  procedure  could  be  iterated  to  produce  a  more  accurate  result,  we 
found  that  iteration  was  not  necessary. 

The  wave  equation  given  by  Eq.  (2-1.29)  may  be  simplified  further  now 
by  comparing  the  relative  sizes  of  the  source  terms.  The  source  terms  involving  the 
pressure  are  ignored  in  this  reasoning  since  the  pressure  is  an  unknown  quantity. 

In  previous  studies  on  thermoacoustic  sources,  the  only  source  term  that 
was  considered  was  the  first  order  term  —  This,  being  the  onlv  first  order 
term,  is  indeed  the  most  significant  driving  term.  It  was  reported  in  a  review 
article  bv  Lvamshev  and  Naugol'nykh2  that  the  term  —  is  the  next  largest 
source  term.  Hence  we  compare  all  of  the  remaining  second  order  terms  to  this 
term,  in  order  to  justify  their  inclusion  or  exclusion  in  our  numerical  analysis  of 
the  wave  equation. 

We  know  that  for  fresh  water  at  room  temperature,  jr  =  3.3  X  10-3. 
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T 

and  —p  -a-  can  be  ignored  when  compared  to  the  one  containing 
wave  equation  can  be  simplified  to: 
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Consider  the  ratio  of  the  second  and  sixth  terms  on  the  right  hand  side 


of  Eq.  (2  2.12): 
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Substituting  Eqs.  (2-2.11)  and  (2-2.8)  for  <$/?  and  7),  we  get 

r>  .  •  ,  1  d(l  f  ,, 


tio  1  =  — ^  , 

2</2  a#  y 


where  the  integral  is  evaluated  over  the  time  interval  during  which  a  section  of 
the  fluid  is  being  heated.  Consider  a  point  in  the  fluid  heated  by  the  laser  beam. 
If  the  beam  is  assumed  to  have  a  circular  cross-section,  then  the  point  is  exposed 
to  heat,  at  most,  for  the  time  it  takes  for  the  laser  beam  to  travel  a  distance 
equal  to  the  beam  diameter.  If  the  diameter  of  the  beam  is  2  cm  and  the  laser  is 
moving  at  Mach  1,  (corresponding  to  our  experimental  setup)  then  the  section  of 
fluid  is  heated  for  a  time  th  =  10-5  s.  Therefore  the  ratio  can  be  written  as 


1  dq  fto+th 
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(2-2.15) 


where  t0  varies  such  that  0  <  <  (tp  -  th).  A  computer  program  has  been 

written  to  calculate  the  maximum  value  of  Ratio  1,  which  was  found  to  be  of 
the  order  of  0.1.  This  justifies  the  elimination  of  the  term  —5 W~qt  >n  comparison 
to  the  term  -jS-Tt  ffrl  .  Note  that  this  conclusion  is  only  valid  if  the  source 
is  moving  at  approximately  Mach  1.  If  the  laser  beam  is  moved  more  slowly, 
then  th  and  Ratio  1  increase.  Since  our  concern  is  with  an  explanation  of  the 
discrepancy  between  theory  and  experiment  observed  for  a  transonic  source,  we 
have  concentrated  on  a  source  moving  at  Mach  1. 

Consider  the  ratio  of  the  third  and  seventh  terms  on  the  right  hand  side 


of  Eq.  (2-2.12): 


Ratio  2 


9/7  rj1  dCp 
1  II  3T  To 


(2-2.16) 


("sing  Eqs.  (2  2.10)  and  (2  2.8).  we  see  that  Ratio  2  =  Ratio  1.  Therefore,  by 


he  same  analvsis  used  above,  the  term  1 '  p‘:y can  be  neglected  when  compared 
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Consider  the  ratio  of  the  seventh  and  sixth  terms  on  the  right  hand  side 


of  Eq.  (2-2.12): 


Ratio  3 


ac P  i 

8o  9T  I  To 

r  33  I 
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which  is  roughly  1.04  X  10-6  for  fresh  water  at  room  temperature.  Thus  it  is  seen 
that  the  seventh  term  can  be  neglected  in  relation  to  the  sixth  one. 

Estimating  the  ratio  of  the  fourth  and  sixth  terms  on  the  right  hand  side 
of  Eq.  (2-2.12)  for  fresh  water, 

8 o  qt6T 


Ratio  4  = 


2  /q  qTt 
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The  order  of  magnitude  of  Ratio  4  justifies  the  exclusion  of  the  fourth  term. 

Consider  the  ratio  of  the  eighth  and  sixth  terms  on  the  right  hand  side 
of  Eq.  (2-2.12): 

Ratio  5  =  — •  (2-2.19) 

Zql  t 

We  may  find  an  expression  for  Ttt  by  differentiating  Eq.  (2-2.8).  Since  6s  = 
(J  qdt)/(p0T0),  this  leads  to  the  conclusion  that  Ratio  5  and  Ratio  1  are  the  same 
which  supports  the  exclusion  of  the  eighth  term  compared  to  the  sixth  one.  Ex¬ 
amining  the  ratios  defined  above,  it  is  obvious  that  the  sixth  term  on  the  right 
hand  side  of  Eq.  (2-2.12)  is  indeed  the  biggest  second  order  term  (excluding  the 
possible  exception  of  the  pressure  dependent  terms,  which  are  evaluated  itera¬ 
tively). 


2-3  NUMERICAL  SOLUTION  OF  THE  WAVE  EQUATION  WITH 
FIRST  AND  SECOND  ORDER  SOURCE  TERMS 

From  the  arguments  presented  above  we  may  rewrite  the  simplified  wave 
equation  as 

V2  _  Pu  _Po(h  _  dfi_  PoPqt  /Wt  f2o.x 

p  cj  C„  C„  dT  n  +  f,0ciCn  +  wlCn  '  1  "  ’  ' 

Computer  programs  to  evaluate  the  pressure-time  response  due  to  each  of  these 
source  terms  independently  and  for  all  of  them  together  were  written.  To  get  the 
total  response,  the  impulse  response  due  to  each  of  the  source  terms  is  computed 
and  added.  This  can  be  done  because  the  wave  equation  is  linear  in  pressure  and 
therefore  superposition  holds. 

Figure  (2-1)  defines  the  geometry  of  the  laser  source-hydrophone  system. 
In  this  figure,  r0,  0o,  and  <f>0  are  the  range  and  angles  defined  by  the  ray  connecting 
the  hydrophone  to  the  first  column  of  water  insonified  by  the  laser.  The  source  is 
assumed  to  be  traversing  across  an  air-water  interface  with  a  speed  of  Mach  1.015. 
Figure  (2-2)  shows  the  pressure  signal  generated  at  a  range  of  4  m  which  is 
attributable  to  the  leading  order  source  term  in  Eq.  (2-3.1).  Figure  (2-3)  shows 
the  pressure  at  the  same  location  due  to  the  second  term  on  the  right  hand  side 
of  Eq.  (2-3.1).  This  pressure  signal  is  attributable  to  the  fact  that  the  coefficient 
of  theimal  expansion  of  water  varies  with  temperature  (and  hence  with  time  for 
our  system).  Note  that  the  pressure  scale  in  Fig.  (2  3)  is  10-3  times  that  in 
Fig.  (2-2).  Figure  (2  4)  shows  the  pressure  signal  at  the  same  location  which 
may  be  attributed  to  the  pressure  dependent  terms  in  Eq.  (2  3.1).  The  scale 
of  the  pressure  axis  is  the  same  as  that  for  Fig.  (2  3).  Note  that  superposition 
of  these  three  plots  produces  a  graph  which  is  identical  to  Fig.  (2  2).  Hence 
we  conclude  that  the  variation  of  thermodynamic  quantities  during  the  sound 
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generating  process  is  not  significant.  This  dispels  all  doubts  about  the  significance 
of  the  second  order  terms  near  Mach  1,  since  the  entire  analysis  is  centered  around 
that  source  speed. 

It  should  be  noted  that  the  analysis  contained  in  this  chapter  addressed 
only  the  effect  of  nonlinearities  in  the  state  equation  of  the  medium.  Linear  ex¬ 
pressions  for  conservation  of  mass  and  momentum  were  employed.  Hence  the 
possibility  of  convective  nonlinearities  contributing  to  the  observed  unusual  be¬ 
havior  of  the  thermoacoustic  source  traveling  at  speeds  near  Mach  1  has  not  been 
eliminated.  However  it  has  been  reported  by  Pierce21  that  the  effect  of  hydro- 
dynamic  nonlinearities  is  negligible. 


CHAPTER  3 


EFFICIENCY  CONSIDERATIONS  FOR  THE 
THERMOACOUSTIC  ENERGY  CONVERSION  PROCESS 


3-1  EFFICIENCY  DEFINITION 

The  question  of  increasing  the  transduction  efficiencies  for  thermoacoustic 
sources  has  been  studied  by  American  and  Soviet  scientists  over  the  past  few  years. 
When  the  concept  of  thermoacoustic  generation  of  sound  was  first  conceived,  the 
experimental  studies  were  limited  to  stationary  sources.  The  main  disadvantage 
of  the  stationary  thermoacoustic  source  is  its  low  energy  conversion  efficiency. 
Figure  (3-1),  taken  from  Muir,  Culbertson,  and  Clynch,5  shows  an  estimate  of 
the  laser  power  needed  to  generate  given  source  levels  at  a  distance  of  1  m.  They 
predicted  that  “megawatts  of  optical  power  would  probably  be  required  in  the 
megahertz  frequency  region,  with  gigawatts  of  power  required  in  the  kilohertz 
region,"  in  order  for  the  transducer  to  be  useful  for  sonar  applications.  This  low 
efficiency  led  to  the  consideration  of  a  moving  thermoacoustic  source. 

It  was  reported  by  Bunkin  et  al.6  in  1978  that  the  amplitude  of  the 
acoustic  pressure  generated  in  the  thermoacoustic  process  could  be  increased  by 
moving  the  laser  beam  at  high  velocity  through  the  fluid.  A  frequency  domain 
approach  was  used  in  their  analysis  which  showed  significantly  increased  sound 
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pressure  levels  for  a  source  moving  at  speeds  near  Mach  1.  Experimental  ver¬ 
ification  of  these  results  was  provided.  However,  in  a  later  article.  Bunkin  et 
al . '  claimed  that  moving  the  laser  beam  does  not  result  in  an  increased  sound 
pressure  level.  Again,  experimental  evidence  was  offered.  Hence,  there  has  been 
great  confusion  over  the  effect  of  moving  a  laser  source. 

In  the  time  domain  approach  developed  by  Bert  helot,  and  Busch- 
Yishniac1*  for  moving  thermoacoustic  sources,  the  acoustic  pressure  is  driven  by 
a  source  that  depends  on  the  time  rate  of  change  of  laser  intensity.  Experiments 
were  conducted  and  the  results  compared  to  analytical  predictions.  It  was  found 
that  there  is  good  agreement  between  theory  and  experiment  except  for  transonic 
motion  of  the  source,  where  an  experimental  enhancement  of  about  25  dB  in  the 
sound  pressure  level  was  observed. 

Research  performed  by  Pierce  and  Ilsieh 18  on  the  optoacoustic  energy 
conversion  efficiency  yielded  an  energy  equation  which  indicated  that  the  rate  of 
acoustic  energy  production  per  unit  volume  is  proportional  to  the  product  of  the 
rate  of  heat  addition  per  unit  volume  and  the  local  acoustic  pressure  existing  at 
the  same  point.  This  led  to  the  conclusion  that  the  laser  energy  deposition  system 
should  be  designed  in  such  a  way  that  a  traveling  acoustic  wave  was  created,  and 
then  its  amplitude  systematically  pumped  up  by  always  adding  energy  at  places 
where  the  acoustic  pressure  has  a  local  maximum.  We  note  that  the  analysis 
presented  in  the  preceding  chapter  has  shown  that  if  one  includes  the  second 
order  terms  associated  with  entropy  generation  in  the  equation  of  state,  then  the 
second  order  source  terms  also  depend  on  the  product  of  pressure  and  another 
t  hermodynamic  variable.  Hence  t  he  pumping  principle  applies  to  t  he  second  order 
source  terms  as  well  as  the  leading  order  term. 

It  was  concluded  in  ( 'hapter  2  t  hat  l  lie  only  term  t  hat  makes  a  significant 
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contribution  to  the  pressure  generated  by  the  thermoacoustic  source  for  our  laser 
system  is  the  leading  order  term.  Therefore  in  the  analysis  presented  in  this 
chapter  we  ignore  all  nonlinear  effects.  With  this  simplification  the  energy  balance 
equation  can  be  written  as 


(5',ou2),  +  ^(e), 


+  V  • (pu) 


To  pq 

poC  po 


(3-1.1) 


The  first  and  second  terms  on  the  left  hand  side  of  Eq.  (3  1.1)  are  the 
kinetic  and  potential  energies  per  unit  time  per  unit  volume  of  the  fluid,  respec¬ 
tively.  The  third  term  is  the  work  done  on  the  fluid  per  unit  time  per  unit  volume. 
The  term  on  the  right  hand  side  of  Eq.  (3-1.1)  is  the  heat  energy  per  unit  time 
per  unit  volume  which  is  converted  to  acoustic  energy.  This  is  in  accordance 
with  the  law  of  conservation  of  energy  which  states  that  the  input  energy  either 
increases  the  stored  energy  or  is  radiated  away.  Clearly,  the  higher  the  fraction  of 
the  available  laser  energy  converted  to  acoustic  energy,  the  more  efficient  the  ther¬ 
moacoustic  process.  Hence  the  efficiency  for  the  process  was  defined  by  Pierce20 
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(3-1.2) 


where  i)  is  the  volume  of  fluid  to  which  heat  is  added  during  the  thermoacoustic 
generation  of  sound. 

The  interesting  property  of  this  efficiency  is  t  hat  t  he  fract  ion  of  converted 
energy  depends  not  only  on  the  heat  energy  adder!  to  the  medium,  but  also  on 
the  acoustic  pressure  of  the  fluid  into  which  it  is  dumped.  If  the  laser  is  inten¬ 
sity  modulated,  as  in  the  experiments  performed  by  Bert  helot  ,12  then  the  heat 
oscillates  in  time,  producing  an  instantaneous  oscillation  of  the  acoustic  pressure. 
From  Eq.  (3  1.2)  it  is  clear  that  this  process  is  not  very  efficient  because  the 
product  of  pressure  and  heat  oscillates  about  a  low  value.  Instead,  it  the  laser  is 
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list'd  to  always  add  heat  to  the  fluid  which  is  already  at  a  pressure  maximum,  then 
the  product  of  pressure  and  heat,  and  therefore  the  energy  conversion  efficiency, 
will  he  high.  If  Pierce's  argument  is  valid,  one  should  see  a  steady  increase  in 
elliciency  with  time'  under  these  conditions. 

Pierce  and  Ilsieh18  have  examined  the  use  of  this  pumping  principle 
in  light  of  currently  available  lasers.  They  have  concluded  that  even  within  the 
limitations  of  existing  laser  technology,  it  should  be  feasible  to  develop  an  airborne 
laser  system  that  will  create  detectable  sound  at  ranges  up  to  10  km  or  greater. 

3-2  NUMERICAL  COMPUTATION  OF  EFFICIENCY  FOR  THE 
THERMOACOUSTIC  PROCESS 

We  have  numerically  evaluated  the  expression  for  efficiency  given  by 
Hq.  (3  1.2)  by  modifying  the  program  written  by  Bert  helot . ,!l  The  linear  first 
order  wave  equation  for  t hermoacoust ic  generation  of  sound,  derived  by  Larson.1' 
is  used  for  calculating  the  pressure  using  a  convolution  approach.  The  heat  input 
per  unit  time  per  unit  volume.  <7,  is  given  by  Fq.  (2  2.5).  The  volume  integral  in 
Kq.  (3  1.2)  is  evaluated  over  the  laser  irradiated  region. 

If  the  radius  of  the  laser  beam  is  a.  the  laser  pulse  duration  rp.  the  speed 
of  the  laser  across  the  water  surface  t\  and  the  optical  coefficient  of  absorption  o . 
then  the  laser  heated  volume  is  approximately  a  parallelepiped  of  dimensions  2 n. 
i'Tp,  and  1/a.  as  shown  in  f  ig.  (3  2).  The  only  irradiated  area  neglected  by  this 
approximation  is  tfiat  correponding  to  two  half-cylinders  extending  outward  from 
the  edges  normal  to  the  //  axis,  (  sing  the  reference  frame  indicated  in  Fig.  (3  2). 
the  lallowing  inequalities  are  valid  for  the  laser  heated  /one: 


Consider  an  element  0[r,y.z)  of  the  fluid  heated  by  the  laser  as  shown 
in  Fig. (3  3).  The  polar  coordinates  for  this  element  are 


=  Jx2  +  y2  +  22  . 


(3-2.1) 


0  =  arccos (z/r)  , 


O  =  arccos 


The  point  O  is  heated  for  a  time  from 


to  tf  — 


x  +  \/a2  —  y 2 


(3  2.2) 
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Therefore  the  time  integral  for  the  efficiency  for  point  O  is  evaluated  with  the 
limits  t ,  to  tj  since  at  any  time  outside  that  range  the  added  heat  (and  thus 
efficiency)  are  zero.  Note  that  these  limits  are  independent  of  2. 

In  order  to  examine  the  efficiency  defined  by  Eq.  (3-1.2)  as  the  source 
moves  at  various  speeds,  we  calculate  p{t)  for  every  point  0  in  the  laser  heated 
zone.  The  product  pq  is  then  accumulated  for  all  such  points  as  a  function  of 
time,  where  q  is  assumed  to  decrease  exponentially  with  depth  2  and  to  have  a 
Gaussian  shading  through  the  beam  cross-section.  This  permits  us  10  calculate 
the  efficiency  as  a  function  of  time  as  the  laser  moves  through  the  volume.  In 
this  manner  we  may  compare  the  efficiencies  for  subsonic  and  transonic  motions 
of  the  laser  beam. 

Note  that  the  Mach  number  of  the  source,  as  seen  by  various  points  in 
the  laser  irradiated  region,  is  not  a  constant.  It  varies  as  M  sin  0,  cos  o,  .  where  M 
is  the  absolute  Mach  number  of  the  source  anil  0,  and  o,  are  the  angles  defining 
t  he  posit  ion  of  a  point  in  t  he  laser  heated  volume  wit  h  respect  to  t  he  ini  I  ial  source 


position.  For  a  constant  value  of  A/,  it  can  be  shown  that  M  sin  0,  cos  d>l  varies 
from  0  to  A/  for  points  in  the  laser  heated  region. 

['lie  computer  program  written  by  Bert  helot 14  for  calculating  the  pressure¬ 
time  response  has  been  mollified  so  as  to  calculate  t  hi'  efficiency  as  a  function  of 
time.  In  this  calculation,  the  system  properties  are  all  fixed,  and  the  efficiency 
is  found  with  the  time  varying  from  rp/ 10  to  r;,  in  steps  of  r;,/10.  The  laser  is 
assumed  to  be  unmodulated  and  identical  to  that  described  in  Chapter  2.  Cal¬ 
culations  of  the  laser-generated  pressure  are  carried  out  at  1000  points  in  the 
laser  heated  volume.  It  is  to  be  noted  that  the  actual  value  of  acoustic  pressure, 
including  its  phase,  is  used  in  the  calculation  of  efficiency. 

figure  (3-4)  shows  the  variation  of  normalized  efficiency  with  time  for 

the  source  moving  in  a  straight  line  at  Mach  0.5.  The  normalized  efficiency  is 

simply  the  efficiency  defined  in  Eq.  (3-1.2)  divided  by  the  constant  30/(puCp0)- 
l  lie  time  axis  corresponds  to  the  time  after  the  laser  lias  been  turned  on  at  which 
the  efficiency  is  calculated.  The  graph  shows  a  rapid  decrease  in  efficiency  with 
time  at  the  beginning  of  the  motion  of  the  laser.  This  is  caused  by  the  fact 
that  the  pressure  signal  has  a  large  negative  portion  immediately  following  the 
positive  peak  (see  Fig. (2  2)).  These  two  peaks  thus  roughly  cancel  one  another 
in  the  integral  over  time  of  the  product  of  /ay.  1  his  large  pressure  dip  is  due  to 
the  presence  of  t  hi'  air-water  interface  which  causes  substantial  diffraction  and 
diffusion.  After  the  high  amplitude  portion  of  the  pressure  signal  is  passed,  there 
is  a  very  modest  gain  in  efficiency. 

figure  (2  o)  is  a  plot  of  normalized  efficiency  versus  time  for  the  source 

moving  near  Mach  1.  Ih  is  plot  indicates  a  trend  similar  to  that  for  Mach  0.5. 

However,  it  is  scon  that  the  maximum  value  of  efficiency  for  the  transonic  source 
is  lower  than  that  for  the  subsonic  source. 
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The  efficiency  as  defined  by  Eq.  (3-1.2)  is  zero  for  a  supersonic  source 
because  the  laser  is  always  heating  a  region  of  water  which  has  not  yet  received 
sound  from  previously  heated  water.  Hence,  we  conclude  from  Figs.  (3  1)  and 
(3  •">)  that  there  is  no  evidence  that  moving  a  laser  source  at  Mach  1  increases 
the  sound  generating  efficiency  compared  to  a  subsonic  source.  In  fact,  there  is 
evidence  to  the  contrary. 


CHAPTER  4 


EXPERIMENTAL  TESTS 


The  time  domain  analysis  presented  by  Berthelot  and  Busch-Vishniac 
inherently  assumes  the  superposition  principle  to  be  valid  in  the  numerical  com¬ 
putation  of  the  acoustic  pressure  generated  by  a  moving  thermoacoustic  source. 
In  their  approach,  the  total  travel  of  the  laser  beam  is  decomposed  in  time. 
Fi  gure  (4-1)  shows  the  thermoacoustic  source  at  various  points  in  time  as  it 
moves  across  the  water  surface  in  a  straight  line.  The  time  interval  between  two 
successive  source  positions  corresponds  to  At,  which  can  be  chosen  at  will  in  or¬ 
der  to  obtain  a  satisfactory  approximation  for  the  continuous  source  motion.  The 
pressure  signal  at  the  receiver  is  computed  by  summing  the  signals  generated  from 
each  of  the  source  positions  during  the  time  the  laser  is  on,  taking  into  account 
the  time  delay  induced  due  to  the  source  motion.  Thus  the  total  pressure  at  the 
receiver  is  viewed  as  a  convolution-type  summation  in  the  time  domain  between 
the  impulse  response  of  the  system  and  the  optoacoustic  source  strength. 

The  method  described  above  was  shown  by  Berthelot  and  Busch- Vishniac 
to  work  well  for  subsonic  and  supersonic  motions  of  the  source.  However,  for  the 
transonic  case,  the  numerical  predictions  seem  to  underestimate  the  measured 
sound  pressure  level  by  about  25  dB.  This  discrepancy  raised  doubts  about  the 
validity  of  the  superposition  principle  for  t  ransonic  motion  of  the  source.  It  was 
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suspected  that  nonlinear  effects  dominate  near  Mach  1.  and  therefore  that  the 
superposition  principle  might  not  he  valid. 

1  his  chapter  attempts  to  clarify  the  questions  raised  about  application 
of  the  superposition  principle.  Specifically,  experimental  data  are  presented  for 
the  pressure- time  response  under  two  conditions:  continuous  motion  of  the 
thermoacoustic  source  for  a  period  of  roughly  1  ms.  and  the  pressure-time 
response  formed  by  superposing  the  signals  generated  by  traversing  nonover¬ 
lapping  portions  of  the  identical  path.  In  the  latter  case,  sufficient  time  is  al¬ 
lowed  to  elapse  between  the  excitation  of  the  path  segments  so  that  interaction  is 
eliminated.  Experiments  were  performed  for  subsonic,  transonic,  and  supersonic 
motions  of  the  thermoacoustic  source. 

4-1  DESCRIPTION  OF  EXPERIMENT  AND  EQUIPMENT  USED 

The  apparatus  used  in  our  study  is  the  same  as  that  used  in  earlier  studies 
on  moving  thermoacoustic  sources  by  scientists  at  Applied  Research  Laboratories. 
I’he  University  of  Texas  at  Austin  (ARL:UT).  The  main  components  include  a 
high  power  laser  source,  intensity  modulation  system,  a  rotating  mirror  with  a 
variable  frequency  power  supply,  a  hydrophone,  and  an  oscilloscope.  Figure  (4-2). 
taken  from  Ref.  22,  is  a  schematic  of  the  experimental  setup.  A  brief  description  of 
the  various  components  is  given  in  he  following  sections.  A  detailed  explanation 
of  the  apparatus  is  given  in  Culbertson.  Chotiros,  and  Berthelot.22 

4-1.1  LASER 

I  lie  laser  used  in  our  experiments  was  procured  from  Apollo  Lasers. 
Inc.  Figure  (1  4)  shows  the  optical  train  configuration  for  this  laser.  The  laser 
cavity  can  accommodate  either  a  Nd:glass  rod  or  a  ruby  rod.  For  simplicity,  a 
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\’d: glass  rod  was  used  throughout  our  experiments.  The  wavelength  of  the  liglit 
emitted  by  the  laser  with  a  glass  rod  is  1.0b  // rn.  Optical  pumping  of  the  glass 
rod  is  achieved  by  a  flash  lamp.  Modulation  of  the  light  intensity  is  affected  by 
polarizing  the  laser  beam  by  a  Brewster  stack  before  it  strikes  a  Pockels  cell  and  a 
prism  analyzer.  Figure  (11),  taken  from  Culbertson,  Chotiros,  and  Bertheiot." 
depicts  the  principle  of  laser  beam  modulation.  The  Pockels  cell  consists  of  an 
electrically  driven  KD*P  crystal  immersed  in  a  protective  dielectric  fluid.  When 
a  voltage1  is  applied  to  the  terminals  of  the  Pockels  cell,  the  plane  of  polarization 
of  light  passing  through  it  is  rotated  by  an  amount  which  is  related  to  the  applied 
voltage.  The  output  of  the  polarizer,  Pockels  cell,  and  analyzer  system  is  given 
by  the  relation 

I  =  Io  sin2  ^7-  ,  (1-1.1) 

“*'1/2 

where  7  is  the  transmitted  light  intensity,  70  is  the  incident  light  intensity,  V  is 
the  applied  voltage,  and  \\/2  is  the  half-wave  voltage  of  the  KD*P  crystal  (6  kY 
for  the  KD*P  crystal  used).  The  normalized  envelope.  70(/),  of  the  laser  intensity 
in  the  experiment  is  described  by  the  relation 

70(0  =  10.8(f/rp)  exp(-5.0//rp)  .  (1-1.2) 

where  rp  is  the  laser  pulse  duration.  Figure  (4-5)  illustrates  the  variation  of  7o(/) 
with  time. 

As  shown  by  Eq.  (4-1.1).  there  is  no  light  output  by  the  modulation  sys¬ 
tem  when  l  =  0.  When  V  =  \  \/i,  the  laser  beam  passes  through  the  modulation 
system  with  its  intensity  unchanged.  With  this  configuration,  then,  the  beam  can 
be  turned  on  by  a  square  pulse  of  amplitude  \  |/2-  The  duration  for  which  the 
beam  remains  in  the  water  is  controlled  by  the  width  of  this  square  pulse. 

I  lie  laser  beam  thus  obtained  passes  through  three  lenses  mounted  on 


the  covered  optical  rail.  The  three  lenses,  two  convex  and  one  concave,  aid  in 
bringing  the  laser  beam  to  a  sharp  focus  on  the  water  surface.  The  diameter  of  the 
beam  and  its  power  were  optimized  to  get  a  sufficiently  narrow  and  powerful  beam 
on  the  water.  A  narrow  beam  minimizes  diffraction  effects  and  a  powerfvd  beam 
generates  a  stronger  pressure  signal.  A  computer  program  was  used  to  calculate 
the  positions  of  the  various  lenses  relative  to  the  glass  rod.  in  order  to  obtain 
acceptable  values  for  the  diameter  and  power  of  the  beam.  The  beam  is  finally 
deflected  onto  the  water  surface  by  a  90°  prism  and  the  rotating  mirror.  A  low 
power  HeNe  laser  is  used  to  align  the  optical  components.  The  rotational  speed 
of  the  hexagonal  mirror  is  adjustable  to  get  linear  source  velocities  corresponding 
to  the  range  Mach  0.1  to  Mach  2.6  relative  to  the  sound  speed  in  water. 

4-1.2  FIRING  CONTROLLER 

The  firing  of  the  laser  is  designed  to  occur  at  a  time  such  that  the  beam  is 
initially  vertical.  The  laser  beam  then  moves  towards  the  receiver  for  the  duration 
of  the  laser  pulse.  This  coordination  is  achieved  by  means  of  a  firing  controller 
designed  raid  built  at  ARL:UT.  Charging  of  the  laser  capacitors  is  controlled 
manually. 

The  main  function  of  the  firing  controller  is  to  generate  a  pair  of  pulses: 
FIRE  and  OUTPUT.  The  FIRE  signal  is  a  pulse  of  approximately  100  fis 
duration  and  has  an  amplitude  of  about  12  V.  This  signal  is  used  to  trigger  the 
laser.  The  angle  at  which  laser  light  appears  on  the  water  can  be  selected  by  a 
pair  of  dipswitches.  The  decimal  numbers  represented  by  the  two  dipswitches  are 
related  to  the  angle.  i.\  made  by  the  beam  with  the  vertical  at  the  beginning  of 
the  laser  sweep,  by  the  following  equation: 


where  the  -f  sign  is  used  when  the  laser  beam  is  moving  away  from  the  laser 
bench.  Similarly,  the  -  sign  is  used  when  the  beam  moves  towards  the  bench. 
Here,  r>  is  a  positive  integer  accounting  for  the  periodicity  of  the  facets  on  the 
mirror  which  is  set  to  zero  for  simplicity.  The  angle  v  can  be  varied  in  steps  of 
120/255  of  a  degree  with  S- bit  dipswitches. 

Fhe  OUTPUT  signal  is  a  standard  TTL  signal  of  about  1  ms  duration, 
which  is  delayed  after  the  FIRE  pulse.  The  length  of  delay  can  be  controlled  bv 
another  dipswitch.  The  OUTPUT  signal  is  used  to  trigger  the  modulating  signal 
produced  by  the  Exact,  signal  generator,  which  forms  the  input  to  the  high  voltage 
burst  generator.  The  delay  between  the  F'IRE  and  OUTPUT  signals  should  be 
at  least  ISO  // s  when  a  Nd:gla.ss  rod  is  used  as  the  lasing  element.  This  is  because 
there  is  an  intrinsic  time  interval  between  the  discharge  of  the  capacitor  bank  and 
the  actual  lasing  of  the  rod.  Fhe  OUTPUT  signal  governs  the  time  at  which  the 
laser  light  appears  on  the  water. 

Fhe  input  to  the  firing  controller  comes  in  the  form  of  a  pair  of  signals 
from  a  photodiode  assembly  placed  inside  the  rotating  mirror  mounting  fixture. 
The  two  signals  are  delayed  in  time,  corresponding  to  a  C)0°  rotation  of  the  mirror. 
Continuous  rotation  of  the  mirror  produces  a  train  of  pulses  with  a  period  equal 
to  the  time  taken  by  the  mirror  to  rotate  by  f>0°.  The  signal  thus  obtained  is 
also  available  as  one  of  the  outputs  from  the  firing  controller.  This  is  labeled 
7  PA  and  is  used  to  accurately  measure  and  monitor  the  speed  of  rotation  of  the 
mirror. 

1-1.3  RECEIVING  AND  ANALYZING  EQUIPMENT 

1  he  pressure  signal  generated  by  the  moving  thermoacoustic  source  is 
received  by  an  11-50  hydrophone  and  viewed  on  a  Nicole!  I()‘)l  digital  oscilloscope. 


I  lie  hydrophone  is  mounted  on  a  vertical  fixture  attached  to  a  mobile  platform  in 
a  fresh  water  tank.  The  horizontal  range  and  depth  of  the  hydrophone,  relative 
to  the  initial  position  of  the  beam  on  the  water,  can  be  selected  at  will  within 
the  constraints  of  the  tank  size.  The  orientation  of  the  hydrophone  is  shown  in 
Fig.  (4-6).  In  this  figure,  .rz  is  the  vertical  plane  and  xy  is  the  horizontal  plane. 
51.  52,  53,  54  represent  the  source  positions  at  the  beginning  of  each  interval 
constituting  the  total  motion  of  the  source.  The  hydrophone.  /?,  is  placed  such 
that  the  beam  moves  towards  it.  Positioning  the  hydrophone  in  this  manner 
ensures  a  reasonably  uniform  hydrophone  response  for  sound  from  all  points  in  the 
laser  heated  region.  The  specifications  of  the  hydrophone  used  in  our  experiments 
are: 

Model  :  I’SRD  type  II  56 

Frequency  range  10  Hz  to  65  kHz 

Free  field  voltage  sensitivity  :  -171  d  13  re  1  V'/// Pa 

Preamplifier  :  transistorized.,  11  dB  gain 

1  he  signal  received  by  the  hydrophone  is  amplified  and  filtered  as  needed 
to  get  a  clean  and  detectable  waveform  on  the  oscilloscope.  A  Burr-Brown  ampli¬ 
fier  with  a  selectable  20-10  dB  gain  is  used  for  amplifying  the  raw  signal,  which 
is  then  filtered  through  a  Kronhite  bandpass  filter. 

I  lie  voltage-time  waveforms  thus  generated  are  viewed  on  the  oscillo¬ 
scope  with  a  sampling  rate  of  500  ns.  The  oscilloscope  is  triggered  by  a  pulse 
labeled  VIUK  S)  AY'  from  the  laser  control  panel.  This  signal  is  generated  at 
the  in.  taut  the  laser  is  fired.  The  waveforms  stored  in  the  oscilloscope's  memory 
are  then  transferred  to  minidisks  for  further  analysis  on  the  mainframe  computer 


T 


4-1.4  CONTROL  PARAMETERS  FOR  THE  EXPERIMENTS 


Ih«’  laser  used  in  our  study  had  a  pulse  duration  of  roughly  1  ms.  The 
maximum  energy  of  the  unmodulated  laser  beam  was  approximately  20  J.  The 


diameter  of  the  beam  on  the  water  surface  was  of  the  order  of  1  cm. 


Experiments  were  performed  for  subsonic,  transonic,  and  supersonic  mo- 
t  ions  of  t  he  source.  The  velocities  of  the  source  were  chosen  to  correspond  to  Mach 
0.7.  1.0.  and  1.3.  The  hydrophone  position  relative  to  that  of  the  source  at  the 
beginning  of  the  laser  sweep  is  shown  in  F'ig.  (4  7).  This  position  was  maintained 
constant  throughout  the  experiments. 

For  a  typical  speed  of  the  source,  the  laser  was  first  scanned  across 
the  water  surface  for  1  ms.  The  resulting  pressure-time  response  was  recorded. 
Next,  the  laser  output  was  tapped  from  four  separate  time  intervals  which,  when 
summed,  constituted  the  laser  output  for  the  case  of  continuous  motion  for  1  ms. 
This  strategy  is  indicated  in  Fig.  (4  8),  which  shows  the  laser  intensity  for  each 
time  interval.  Each  interval  is  roughly  250  /ts  in  duration.  The  three  dipswitches 
in  the  firing  controller  are  set  to  the  correct  delay  between  the  firing  of  the  laser 
ind  the  first  appearance  of  laser  light  on  the  water;  the  angle  of  the  beam  is 


set  to  be  vertical  when  it  first  strikes  the  water.  A  total  of  four  waveforms  is 


detained  for  any  Mach  number  of  the  source.  The  dclav  between  measurements 


of  these  waveforms  ensures  that  the  signals  they  generate  are  independent.  The 
lour  waveforms,  constituting  the  complete  sweep  of  the  beam,  are  added  together 
and  compared  to  the  waveform  obtained  by  a  single  laser  sweep  of  1  ms  duration. 
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1-2  EXPERIMENTAL  RESULTS 

Figures  (1  9)  ami  (4  10)  show  the  sum  of  the  four  pieces  of  the  pressure 
signal  and  the  pressure  signal  generated  by  a  single  sweep  of  the  transonic  laser 
beam,  respectively.  I  he  signal  detected  by  the  hydrophone  was  captured  on  disk 
without  any  amplification  or  filtering.  Note  that  there  is  good  agreement  between 
the  two  pressure  signals.  The  slight  variation  in  amplitude  may  be  attributed  to 
the  presence  of  noise  in  the  vicinity  of  the  hydrophone.  This  noise  accumulates 
in  the  sum  waveform,  thereby  increasing  its  amplitude. 

Figure  (4  11)  shows  the  composite  of  the  four  separate  waveforms  con¬ 
stituting  the  entire  motion  of  the  source  moving  at  Mach  0.7.  Figure  (4-12) 
illustrates  the  corresponding  pressure-time  response  for  the  single  sweep  of  the 
beam  for  1  ms.  The  signals  obtained  for  this  case  were  amplified  with  a  gain 
of  40  dB  and  then  filtered,  the  band  passed  being  11.5  kHz  to  100  kHz.  Clearly. 
Fig.  (4-11)  indicates  a  higher  pressure  amplitude  compared  to  Fig.  (4-12).  Tins 
discrepancy  might  be  due  in  part  to  the  significant  contribution  to  the  pressure 
signal  by  the  amplified  noise  when  the  four  waveforms  are  added. 

Similar  results  are  noticed  for  the  supersonic  motion  of  the  source,  but 
the  agreement  between  the  composite  shown  in  Fig.  (4-13)  and  the  single  sweep 
pressure  signal  in  Fig.  (4-14)  for  Mach  1.3  is  better  than  that  for  Mach  0.7. 

I  he  gain  and  the  filter  band  for  the  supersonic  case  are  identical  to  those  for  the 
subsonic  case. 

A  quantitative  comparison  of  the  experimental  data  can  be  obtained 
by  calculating  the  rms  percent  deviation  from  superposition  for  various  Mach 
numbers  of  the  source.  It  is  found  that  the  rms  deviation  is  about  207%  for 
Mach  0.7.  255'/  for  Mach  1.3,  and  59%  for  Mach  1.0. 


The  experimental  data  presented  generally  support  the  validity  of  the 
superposition  principle  for  subsonic,  transonic,  and  supersonic  motions  of  the 
thermoacoustic  source  on  the  water  surface.  In  particular,  we  see  very  good  agree¬ 
ment  between  the  composite  signal  and  the  continuous  signal  when  the  source  is 
moving  at  Mach  1,  as  indicated  by  the  rms  percent  deviation.  From  this  we  con¬ 
clude  that  there  is  no  evidence  to  support  the  supposition  that  nonlinear  effects 
are  the  cause  of  the  amplitude  magnification  observed  for  a  thermoacoustic  source 
traveling  at  Mach  1. 


CHAPTER  5 


SUMMARY  AND  CONCLUSIONS 


The  main  objective  of  this  part  of  the  study  was  to  investigate  the  pre¬ 
vailing  theories  concerning  the  causes  of  increased  sound  generation  for  a  transon- 
ically  moving  thermoacoustic  source.  The  approach  employed  three  techniques. 
First  a  wave  equation  was  derived  in  which  the  leading  nonlinearities  were  re¬ 
tained  in  the  equation  of  state.  Second  the  efficiency  defined  by  Pierce  and  Hsieh 
was  numerically  determined.  Third  an  experimental  check  of  the  validity  of  the 
superposition  principle  was  conducted. 

The  wave  equation  for  sound  generation  by  the  thermal  mechanism  was 
rederived  including  the  most  important  second  order  term  in  the  equation  of 
state.  The  effects  of  time  varying  the  thermodynamic  properties  of  the  medium 
were  thus  considered.  The  new  forcing  terms  in  the  wave  equation  include  the 
effects  of  a  time  varying  specific  heat,  coefficient  of  thermal  expansion,  acoustic 
pressure,  and  temperature.  The  wave  equation  was  solved  numerically  with  the 
linear  source  term  and  the  second  order  terms  involving  the  coefficient  of  thermal 
expansion  and  the  acoustic  pressure.  The  ensuing  analysis  of  the  sound  field  gen¬ 
erated  by  the  source  does  not  indicate  a  significant  variation  over  that  obtained 
by  earlier  studies  on  transonic  thermoacoustic  sources.  Therefore  it  can  be  con¬ 
cluded  that  the  variation  with  time  of  the  pertinent  thermodynamic  quantities. 
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(  p  and  d,  has  little  effect  on  the  pressure  generated  by  a  moving  thermoacoustic 
source.  The  linear  wave  equation  derived  earlier  by  Larson13  with  a  single  first 
order  driving  term  is  perfectly  valid  for  studies  on  thermoacoustic  sources  over 
the  entire  range  of  source  velocities. 

The  concept  of  efficiency  for  moving  thermoacoustic  sources  has  been 
examined.  The  efficiency  of  the  thermoacoustic  process  was  defined  by  Pierce 
as  the  fraction  of  the  heat  energy  that  is  converted  to  acoustic  energy.  It  was 
suggested  that  moving  the  source  at  Mach  1  would  significantly  improve  the  trans¬ 
duction  efficiency  of  the  thermoacoustic  process  because  of  the  resulting  dumping 
of  heat  in  regions  of  pressure  maxima.  As  a  direct  consequence  of  this  princi¬ 
ple,  we  should  expect  a  monotonically  increasing  value  of  efficiency  with  time 
during  the  interval  of  heat  addition.  The  integrals  constituting  the  definition  of 
efficiency  were  evaluated  numerically  and  the  resulting  value  of  efficiency  plotted 
as  a  function  of  time.  These  plots  indicate  a  rapidly  decreasing  efficiency  with 
time  at  the  beginning  of  the  signal,  for  any  Mach  number  of  the  source.  This  is 
because  the  effect  of  the  air-water  interface  is  to  diffuse  and  diffract  the  sound 
wave,  thus  tending  to  decrease  the  efficiency  of  the  process.  For  later  portions  of 
the  signal  there  is  a  modest  increase  in  efficiency.  The  plots  also  show  that  the 
efficiency  is  higher  for  subsonic  motion  of  the  source  than  for  transonic  motion. 

The  experimental  studies  performed  in  this  research  support  the  assump¬ 
tion  that  the  superposition  principle  applies  for  moving  thermoacoustic  sources. 
The  superposition  principle  has  been  verified  for  subsonic,  transonic,  and  su¬ 
personic  motions  of  the  source.  This  supports  the  theory  that  the  linear  wave 
equation  is  applicable  to  this  process.  Tims  we  see  no  evidence  to  support  the  pre¬ 
diction  that  nonlinear  effects  are  important  for  the  transonically  moving  source. 

Reviewing  the  arguments  presented  in  this  part  of  the  report,  it  is  clear 
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that  none  of  the  theories  considered  was  able  to  explain  the  observed  enhancement 
of  the  pressure  signal  when  the  source  moved  at  Mach  1.  This  is  because  all  of 
them  concentrate  on  the  significance  of  nonlinear  effects.  The  conclusions  reached 
in  this  analysis  indicate  the  validity  of  the  linear  theory. 

The  results  of  this  investigation,  in  conjunction  with  the  articles  pub¬ 
lished  in  Soviet  journals,  suggest  that  a  more  accurate  model  of  the  laser  intensity 
distribution  across  the  beam  cross-section  could  provide  clues  to  the  unpredictable 
behavior  near  Mach  1.  Presently,  it  is  assumed  that  the  intensity  has  a  (laussian 
shading  through  the  beam  cross-section. 


CHAPTER  6 


INTRODUCTION 


The  work  on  thermoacoustic  sources  which  has  been  reported  in  the  liter¬ 
ature  in  the  last  decade  has  typically  concentrated  on  the  efficiency  and  directivity 
of  the  sound  generated  in  this  process.  Although  important,  these  measures  ob¬ 
scure  the  fact  that  the  sound  signal  generated  by  a  moving  source  varies  in  form 
as  the  source-receiver  geometry  changes.  In  the  research  reported  here,  we  study 
the  acoustic  signal  generated  by  the  moving  laser  in  detail.  Properties  of  the  sig¬ 
nal  are  identified  and  the  manner  in  which  these  properties  change  with  respect 
to  source  and  receiver  parameters  is  quantified.  This  information  identifies  those 
aspects  of  the  sound  signal  that  are  most  and  least  susceptible  to  change.  If  the 
thermoacoustic  source  is  used  for  underwater  communication  in  the  future,  the 
statistics  compiled  here  should  be  quite  useful. 

In  the  work  described  here  a  FORTRAN  computer  program  called  MTS. 
written  by  Yves  II.  Berthelot,19  is  used  to  numerically  predict  the  received  acoustic 
signal.  This  computer  program  is  based  on  the  time  domain  approach  discussed 
previously.  Various  properties  of  the  received  acoustic  signal,  such  as  maximum 
and  minimum  periods  and  the  duration  of  the  received  signal,  are  examined  as 
a  function  of  the  source  parameters,  such  as  velocity  and  source- receiver  range 
and  angle.  This  enables  us  to  quantitatively  determine  which  signal  properties 
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SYSTEM  PARAMETERS 


The  laser  assumed  to  be  generating  the  sound  in  this  part  of  the  study 
is  a  Neodymiurmglass  laser  identical  to  the  laser  discussed  earlier.  This  laser  has 
a  wavelength  of  1.06  /mi,  and  a  corresponding  optical  coefficient  of  absorption,  a. 
of  13.7  Np/m  in  fresh  water  where  the  speed  of  sound,  c,  is  1486  m/s.  The  laser 
pulse  duration  was  set  to  0.8  ms,  while  the  laser  beam  radius  was  set  to  0.5  cm. 
A  Gaussian  intensity  distribution  across  the  laser  beam  was  assumed. 

The  parameters  of  the  source  which  are  varied  may  be  classed  into  two 
categories:  those  associated  with  the  source,  and  those  associated  with  the  source- 
receiver  geometry.  The  parameters  of  the  source  which  are  varied  are  as  follows: 
(1)  the  modulation  frequency,  /,  of  the  laser  intensity,  shown  in  Fig.  7-1,  and  (2) 
the  velocity,  v,  or  Mach  number,  M  =  v/c ,  of  the  laser  beam  on  the  surface  of 
the  water,  as  shown  in  Fig.  7-2. 


The  parameters  associated  with  the  source- receiver  geometry  which  are 


varied  are  as  follows:  (1)  the  initial  distance  between  the  source  and  the  receiver. 


r0.  (2)  the  initial  horizontal  angle.  <p,.  and  (3)  the  depth  of  the  receiver,  h.  All  of 


these  parameters  are  shown  in  Fig.  7  2. 

Table  7  1  shows  the  values  of  the  parameters  which  were  used  in  the 


study.  All  possible  combinations  of  the  parameters  were  used,  resulting  in  a  total 
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Figure  7  1  Monci. a i kf>  Laser  Intensity 
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of  675  cases.  These  values  of  the  parameters  were  carefully  chosen  to  cover  the 
whole  range  of  physically  realizable  situations.  The  only  limitation  on  them  is 
that  the  hydrophone  depth,  h,  is  limited  to  a  maximum  value  of  approximately 
54  m,  since  it  appears  in  the  argument  of  the  hyperbolic  sine  function  sinhfo  *  h) 
in  the  program  MTS,  and  the  argument  itself  is  limited  to  a  maximum  value  of 
approximately  742  on  our  computer.  However,  this  limitation  is  a  minor  one. 
since  depths  greater  than  54  m  are  unlikely  to  be  used  for  hydrophones. 

7-1  EXAMINED  PROPERTIES  OF  THE  RECEIVED  SIGNAL 

The  examined  properties  of  the  received  acoustic  signal  may  be  placed 
into  two  categories:  time  related  properties  and  pressure  related  properties. 

The  time  related  properties  of  the  acoustic  signal  which  were  numerically 
investigated  are  as  follows:!  1)  the  duration  of  the  received  signal,  TD.  (2)  the 
time  at  which  the  absolute  peak  pressure  occurs,  TPK.  with  zero  time  defined 
as  the  time  at  which  the  received  signal  begins,  (3)  the  maximum  period  in  the 
signal,  TM AX,  found  from  zero  crossings  of  the  signal,  (4)  the  minimum  period. 
I'M  IN,  (5)  the  average  of  all  periods,  TAX  E,  and  (G)  the  ratio  of  the  time 
at  which  the  peak  pressure  occurs  to  the  duration  time,  TPh/TD.  As  will  be 
discussed,  T PK /TD  is  a  measure  of  the  amount  of  time  inversion  of  the  signal. 

These  properties  are  shown  in  Figs.  7  3(a)  and  (b).  The  beginning  and 
end  of  the  received  signal  are  defined  as  the  points  where  the  sound  pressure 
level  of  the  signal.  Lp ,  has  dropped  to  40  dll  below  the  of  the  absolute  peak 
pressure.  Because  of  noise  inherent  in  a  real  signal,  computing  the  first  and  last 
periods  would  be  difficult.  Therefore  they  were  not  considered  in  this  stud\  mdess 
the  signal  contained  only  three  periods.  In  this  study  that  only  occurs  when  the 
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Mach  number  as  seen  by  the  receiver.  RM,  is  approximately  equal  to  one.  The 
Mach  number  observed  by  the  hydrophone  is  defined  by  RM  =  M  sin  0^  cos  d>0 
where  0o  —  arccos(/;/r,)  (Fig.  7-2  shows  how  6q  and  <po  are  defined).  Note  that 
RM  only  equals  the  Mach  number  of  the  laser  beam  when  the  source  is  moved 
in  a  straight  line  directly  toward  the  hydrophone. 

The  pressure  related  properties  which  we  examined  are  as  follows:  (1)  the 
absolute  peak  pressure,  P PI\  ,  as  shown  in  Fig.  7-3(b),  (2)  the  root  mean  square 
pressure.  P RMS.  of  the  acoustic  signal,  and  (3)  the  pressure  ratio,  PPK/PRMS. 

In  order  to  obtain  the  pressure  and  time  properties  just  described,  approx¬ 
imately  110  FORTRAN  lines  were  added  to  the  already  existing  MTS  program. 
This  code  solves  for  the  pressure  signal  as  a  function  of  time  using  the  previously 
described  convolution-like  approach.  In  general,  with  a  time  step  of  0.671  /*s, 
each  case  took  approximately  16  minutes  to  run  on  a  CYBER  830  Control  Data 
Corporation  computer  system  located  at  ARL:UT. 


CHAPTER  8 


RESULTS  OF  NUMERICAL  STUDY 


In  this  chapter  we  present  the  numerically  obtained  results  for  the  vari¬ 
ation  of  the  signal  properties  as  the  source  and  geometric  properties  of  the  laser 
sound  system  are  varied.  Each  signal  parameter  is  discussed  independently  in 
terms  of  its  sensitivity  to  changes  in  the  system  characteristics. 

8-1  TIME  RELATED  PROPERTIES 

8-1.1  DURATION  TIME,  TD 

The  numerical  results  obtained  for  TD  are  shown  in  Fig.  8-1.  Figure  8-1 
is  typical  of  the  figures  that  will  be  presented  in  this  chapter.  A  total  of  seven 
graphs  are  presented,  each  displaying  TD  versus  bo-  In  each  graph  there  are  five 
curves,  each  corresponding  to  a  different  source  velocity. 

Figures  8-l(a)-(g)  show  that  TD  is  quite  sensitive  to  changes  in  the 
initial  horizontal  angle,  bo,  and  the  Mach  number,  A/,  of  the  laser  beam  on  the 
surface  of  the  water.  Both  of  these  variables  affect  the  apparent  Mach  number  of 
t  he  source  as  observed  by  the  hydrophone,  so  the  sensitivity  of  TD  to  them  comes 
as  no  surprise.  However,  note  that  TD  is  not  sensitive  to  changes  in  the  angle 
bo  when  Xt  =  0  since  the  stationary  source  radiates  sound  omnidirectionally.  TD 


generally  becomes  more  sensitive  to  changes  in  (f>0  as  A/  increases  from  0  to  2 
since  the  source  tends  to  become  more  directive  as  the  source  speed  increases. 

We  also  note  that  TD  is  least  sensitive  to  changes  in  M  when  d>o  =  90°  so 
that  the  laser  beam  moves  in  a  path  normal  to  the  line  connecting  the  hydrophone 
and  the  starting  point  of  the  laser.  The  apparent  Mach  number,  RM ,  under  these 
conditions,  begins  at  zero  and  slowly  increases  as  <f>,  changes.  From  the  data  for 
the  cases  with  4>q  =  90°  we  see  that  TD  increases  slightly  as  M  increases  from 
0  to  2,  since  the  range  of  angles  of  <f>i  increases,  but  the  difference  between  TD 
for  M  =  0  and  M  =  2  is  only  on  the  order  of  2%.  Note  that  T D  generally  becomes 
more  sensitive  to  changes  in  M  as  <f>o  decreases  from  90°  to  0°,  and  increases  from 
90°  to  180°.  Restated,  TD  is  most  sensitive  to  changes  in  M  when  the  laser  is 
moving  directly  towards  the  receiver,  or  directly  away  from  the  receiver. 

Viewing  Fig.  8-1,  we  conclude  that  the  d  .ration  of  the  signal  increases 
monotonically  as  the  laser  moves  away  from  the  receiver,  that  is,  when  <f>,  is  in  the 
range  of  90-180°,  and  that  the  maximum  duration  of  the  received  signal  occurs 
when  M  =  2  and  d>0  =  180°.  This  result  can  be  explained'by  the  fact  that  if  the 
source  is  moving  away  from  the  hydrophone,  then  the  distance  which  the  last 
wavelet  must  travel  from  source  to  receiver  increases  as  the  velocity  of  the  source 
increases,  and  is  greatest  when  the  laser  is  moving  directly  away  from  the  receiver. 

Finally,  an  explanation  can  be  given  of  the  manner  in  which  TD  varies 
as  the  laser  moves  towards  the  receiver,  that  is,  when  <f> 0  is  in  the  range  of  0 
90°.  In  this  region,  the  received  signal  becomes  more  compressed  in  time  as  RM 
increases  from  0  to  1.  and  will  then  become  less  compressed  as  RM  increases 
further.  Pile  shortest  duration  of  the  signal  should  then  occur  near  /?.!/  =  1.  since 
for  this  case  the  whole  signal  arrives  at  the  receiver  essentially  at  the  same  time. 
Ihe  duration  of  the  signal  increases  as  RM  becomes  greater  than  one,  because 


the  distance  which  the  last  wavelet  must  travel  to  reach  the  receiver  decreases  as 
the  velocity  of  the  source  increases.  Therefore,  the  last  wavelet  emitted  by  the 
source  arrives  at  the  receiver  increasingly  earlier  than  the  first  wavelet  emitted 
by  the  source  as  the  velocity  of  the  source  increases.  'The  result  is  that  the  signal 
duration  increases  with  the  source  velocity.  Proof  of  this  is  shown  in  Fig.  8-2. 
which  displays  I  D  versus  RM  for  one  set  of  system  parameters. 

Figures  8  1(a)  (c)  represent  identical  cases  except  for  the  depth  of  the 
receiver.  From  these  three  graphs  one  sees  that  TD  is  very  insensitive  to  changes 
in  the  depth,  /;,  of  the  receiver.  The  data  for  these  cases  reveal  that  I'D  increases 
slightly  as  h  increases.  The  increase  between  the  cases  for  which  h  =  5  m  and 
h  =  -10  m  is  generally  on  the  order  of  0.3%. 

The  cases  used  to  generate  Figs.  8-l(e)-(g)  are  identical  except  for  the 
modulation  frequency  of  the  laser,  /.  prom  these  figures  we  conclude  that  Tl)  is 
very  insensitive  to  changes  in  the  modulation  frequency.  The  data  for  these  cases 
show  that  the  duration  of  the  received  signal  increases  slightly  as  the  modulation 
frequency  increases.  The  increase  between  the  cases  for  which  /  =  r)  kHz  ami 
/  =  00  kHz  was  on  the  order  of  0.6%.  This  slight  increase  in  TD  as  /  increases 
car.  be  explained  as  follows:  Bert  helot  and  Busrh-Y ishniac 12  have  shown  that  the 
pressure  amplitude  of  the  received  signal  increases  as  /  increases  (as  will  also  be 
shown  later  by  the  pressure  related  properties).  Therefore  the  beginning  and  end 
of  t  lie  received  signal  are  detec  ted  eai  her  and  later,  respec  t  i\elv.  by  t  lie  r<-<  fiver, 
causing  an  increase  in  /  D. 

Figures  S  1(a).  (d).  and  (e)  were  generated  horn  identical  sjtu.itjons 
except  lor  the  source  to  receiver  initial  range-.  I  h<-\  -dime  that  I  l>  o  al-o  \ci\ 
insensitive  to  changes  in  the  initial  distance  between  the  vuiur  ami  receive)  .  rM 


I  he-  rest  1 1 1  s  lor  t  hese  rases  show  I  hat  t  he  <  1  u  ra  t  it  mi  o(  t  lie  i  ee  ei  \e<  |  signal  <  !<•<  1 1  a 
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slightly  as  r0  increases.  The  difference  between  the  cases  for  which  r0  =  100  m  and 
r0  =  2000  m  is  generally  on  the  order  of  3 %.  This  result  can  be  explained  by  the 
fact  that  spherical  spreading  will  decrease  the  pressure  amplitude  of  the  received 
signal,  and  that  spherical  spreading  increases  as  r0  increases.  Therefore  the  begin¬ 
ning  and  end  of  the  received  signal  are  detected  later  and  earlier,  respectively,  by 
the  receiver,  and  thus  the  decrease  in  TD. 

From  these  results  we  conclude  that  the  duration  of  the  sound  signal 
generated  by  the  thermoacoustic  source  is  insensitive  to  changes  in  r0 ,  /,  and  h, 
but  is  dependent  on  the  apparent  speed  of  the  source  as  seen  by  the  receiver.  In 
terms  of  geometrical  parameters,  this  means  that  the  duration  of  the  signal  is 
highly  dependent  on  <jf>o  and  A/. 

8-1.2  TIME  OF  PEAK  PRESSURE,  TPK 

Figure  8-3  shows  the  results  for  TP1\ ,  the  time  at  which  the  peak  pres¬ 
sure  is  observed.  Like  the  time  property  TD,  TPK  is  quite  sensitive  to  changes 
in  00  and  A/,  as  shown  by  Figs.  8-3(a)-(g).  For  a  stationary  source  the  graphs  in 
Fig.  8-3  show  that  TPK ,  like  TD,  is  not  a  function  of  d>o  since  the  same  signal 
would  be  observed  at  all  angles.  However,  unlike  TD.  TPK  does  not  increase 
monotonically  as  the  laser  moves  away  from  the  receiver,  although  there  is  a  pat¬ 
tern  in  the  way  TPK  varies  with  A/  when  ct>()  =  0° .  Figure  8  4.  which  displays 
TPK  versus  M  for  a  specific  case,  more  clearly  reveals  this  pattern.  It  is  im¬ 
portant  to  note  that  for  this  specific  case  PM  -  A/  sin  (h>  cos  <p0  is  approximately 
equal  to  A/,  since  <■>„  =  0°  and  0{)  m  90°.  Figure  S  1  shows  that  /'PA  decreases 
as  A/  goes  from  0  to  1.  and  increases  as  A/  goes  from  1  to  2.  Ibis  is  a  result  of 
the  signal  being  compressed  in  I  ime  ( Doppler  shitted )  when  A  /  =  I  is  approached 
from  below  and  above.  This  effect  of  time  straining  would  lead  one  to  expect 
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Figure  8  4  Time  of  Peak  Pressure  versus  Mach  Number 
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Fig.  8-3  to  be  roughly  symmetric  about  A/  =  l.  However,  a  second  effect  result¬ 
ing  from  inversion  of  the  signal  when  M  >  1  skews  the  curve.  In  our  simulated 
signals  from  a  stationary  source,  the  peak  pressure  occurs  relatively  early.  If  the 
source  is  moved  rapidly,  so  that  the  entire  source  line  appears  to  the  receiver  to 
be  moving  with  M  >  1,  then  the  signal  is  inverted.  This  has  the  effect  of  delaying 
the  time  after  signal  onset  at  which  the  peak  pressure  is  observed.  From  Fig.  8-4 
we  conclude  that  the  time  inversion  is  more  important  above  M  =  1  than  the  time 
straining. 

The  cases  used  to  generate  Figs.  8-3(a)-(c)  are  identical  except  for  the 
depth  of  the  receiver.  On  comparing  these  figures  it  can  be  concluded  that  TPK 
is  sensitive  to  changes  in  h,  and  more  so  when  the  velocity  of  the  laser  on  the 
surface  of  the  water,  t>,  is  greater  than  the  speed  of  sound  in  water,  c.  However, 
note  that  TPK  is  not  sensitive  to  changes  in  h  when  A/  =  0. 

Figures  8-3(a),  (d),  and  (e)  were  generated  from  identical  situations 
except  for  the  source-to-receiver  range.  Except  when  M  =0,  TPK  is  sensitive  to 
changes  in  r0,  and  more  so  when  v  >  c.  Note  that  varying  either  r0  or  h  has  the 
effect  of  changing  the  apparent  source  speed.  Hence  it  is  not  surprising  that  the 
pattern  of  signal  variations  is  similar  in  the  two  cases. 

Figures  8-3(e)-(g)  represent  identical  cases  except  for  the  modulation 
frequency  of  the  laser.  They  show  that  TPK  is  sensitive  to  changes  in  f.  even 
when  A/ =  0.  The  fact  that  the  frequency  of  modulation  affects  TPK  even  for 
a  stationary  source  is  primarily  a  low  frequency  effect.  At  very  high  modula¬ 
tion  frequencies,  one  would  expect  that  the  laser  intensity  envelope  peak  would 
determine  the  sound  peak  time  regardless  of  the  modulation. 

From  these  results  we  coin  lude  that  the  time  of  the  peak  pressure  is 
sensitive  to  changes  in  g>().  M,  h.  f.  and  r„. 
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8-1.3  TIME  INVERSION  PROPERTY,  TPK/TD 


For  our  laser,  the  peak  amplitude  of  the  laser  intensity  occurs  relatively 
early.  Therefore,  the  ratio  TPK/TD  can  be  used  as  a  measure  of  the  degree  of 
time  inversion  of  the  received  acoustic  signal.  However,  it  is  important  to  note 
that  for  different  laser  intensity  profiles  TPK/TD  may  not  he  a  measure  of  time 
inversion.  For  example,  in  a  symmetric  pulse  with  the  peak  occurring  exactly 
ft.  the  midpoint  of  the  pulse,  TPK/TD  would  remain  constant  at  a  value  of  0.1 
regardless  of  the  degree  of  time  inversion. 

Since  time  inversion  of  the  signal  can  only  occur  if  the  source  is  moving, 
no  time  inversion  takes  place  if  M  =  0.  Therefore  the  ratio  of  I'PK/'l'D  for 
M  =  0  can  he  used  as  a  reference  for  measuring  the  level  of  time  inversion  of 
the  received  signal.  Figure  8-5  shows  the  results  for  TPK/TD.  My  comparing 
Figs.  8  5(a)  (g)  one  sees  that  TPK/TD  is  generally  not  sensitive  to  changes  in 
On  when  ;•  <  r,  hut  is  quite  sensitive  to  change's  in  on  when  r  >c  arnl  the  laser  is 
moving  towards  the  receiver.  Also,  note*  that  TPK/I  D  is  gene-rail)  not  se-tiMtive 
to  changes  in  \f  when  the*  laser  is  moving  away  from  the  reveive-r.  hut  is  sensitive 
to  changes  in  M  when  the-  laser  is  moving  towards  the  reveive-r.  ('specially  as  M 
goes  from  0.5  to  1.5.  In  other  words,  the  received  signal  is  heavily  time'  inverted 
whe'ii  the'  lase'r  velocity  is  greate'r  than  the-  sound  spewed,  and  is  moving  toward- 
the  re  vei  ver.  Furt  !ie*niie>re*.  the  level  uf  time  inve'rsjon  F  highlv  dependent  on  y, . 
and  M 
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is  fairly  insensitive  to  changes  in  /.  As  in  the  case  of  TPK,  any  variation  in 
TPK/TD  would  be  expected  to  be  a  low  frequency  affect  associated  with  delays 
in  the  peak  of  the  laser  intensity  due  to  modulation. 

Figures  8-5(a),  (d).  and  (e)  were  generated  from  identical  situations 
except  for  the  source-to-receiver  range.  They  sfv  w  that  TPK/TD  is  also  fairly 
insensitive  to  changes  in  r0. 

From  these  results  we  conclude  that  the  time  property  TPK/TD  is 
primarily  sensitive  to  changes  in  <f> o  and  A/,  so  long  as  M  >  1.  Therefore  the 
level  of  time  inversion  of  the  signal  is  very  sensitive  to  changes  in  4>0  and  M. 

8-1.4  MINIMUM  PERIOD,  T M IN 

Figure  8-6  shows  the  results  for  TMIN,  the  minimum  period  in  the 
received  signal.  Figures  8 -6(a)- (g)  show  that  TMIN  is  quite  sensitive  to  changes 
in  <t>0  and  M.  However,  notice  that  TMIN  is  not  sensitive  to  changes  in  <p0  when 
the  laser  is  stationary.  Also,  notice  that  when  </>o  =  90°,  TMIN  is  the  same 
regardless  of  the  velocity  of  the  laser.  The  latter  result  stems  from  the  fact  that, 
when  <£o  =  90°,  the  Mach  number  as  seen  by  the  receiver  is  initially  equal  to  zero 
and  only  increases  slightly.  Therefore  the  signal  almost  appears  to  be  generated 
by  a  stationary  source.  Nevertheless,  one  w'ould  expect  T MIN  to  increase  slightly 
as  M  increases  from  0  to  2  at  0O  =  9 0°,  since  the  signal  will  be  slightly  stretched 
in  time  when  the  laser  is  moving  away  from  the  receiver;  the  data  in  which  results 
were  printed  up  to  the  fourth  significant  digit  do  not  show  this.  It  is  most  likely 
that  the  difference  was  simply  too  small  to  detect. 

Figures  8-6(a)-(g)  show  that  TMIN  is  particularly  sensitive  to  changes 
in  the  velocity  of  the  laser  when  it  is  moving  directly  towards  and  away  from 


show  that  TMIN  increases  monotonically  (with  respect  to  M)  when  the  laser  is 
moving  away  from  the  receiver.  The  largest  value  of  T M IN  in  all  cases  occurs 
when  <po  =  180°  and  M  =  2.  As  discussed  earlier,  this  case  produces  the  greatest 
stretching  of  the  received  signal,  and  therefore  the  largest  TMIN. 

Figures  8-6(a)-(c)  were  generated  from  identical  cases  except  for  the 
depth  of  the  receiver.  They  show  that  TMIN  is  sensitive  to  changes  in  h,  partic¬ 
ularly  when  the  source  is  moving  towards  the  receiver,  and  M  equals  1  or  1.5. 

The  cases  used  to  generate  Figs.  8-6(e)-(f)  are  identical  except  for  the 
modulation  frequency  of  the  laser.  As  would  be  expected,  these  figures  show 
that  TMIN  decreases  as  the  frequency  increases.  Also,  the  data  for  these  cases 
indicate  that  TMIN  acts  almost  as  if  it  were  inversely  proportional  to  /,  regard¬ 


less  of  the  Mach  number. 


Figures  8-6(a),  (d),  and  (e)  represent  identical  situations  except  for  the 
source- to- receiver  range.  These  figures  show  that  TMIN  is  not  particularly  sensi¬ 
tive  to  changes  in  r0. 

From  these  results  we  conclude  that  the  time  property  TMIN  is  sensitive 
to  changes  in  /?,  /,  M,  and  <An  In  addition,  the  manner  in  which  TMIN  varies 
with  respect  to  /,  M,  and  <£0  is  quite  evident. 


8-1.5  MAXIMUM  PERIOD,  T MAX 


Figure  8-7  shows  the  results  for  TM AX,  the  maximum  period  in  the 
received  signal.  Figures  8-7(a)-(g)  show  that  TM AX  is  very  sensitive  to  changes 
in  <f>  o  and  M.  As  for  the  previous  properties,  TM  AX  is  not  sensitive  to  changes 
in  0o  when  the  laser  is  stationary.  All  the  graphs  in  Fig.  8-7  indicate  that  TM  AX 
is  least  sensitive  to  changes  in  M  when  <^>o  =  90° .  Again,  for  this  case  the  Mach 
number  as  seen  by  the  receiver  is  nearly  equal  to  zero  for  the  durat  ion  of  the  pulse: 


therefore  the  sound  appears  to  be  generated  by  an  almost  stationary  source.  How¬ 
ever,  the  data  for  these  cases  do  show  that  TMAX  generally  increases  slightly 
as  the  velocity  of  the  laser  increases.  Like  TMIN,  TMAX  increases  monotoni 
tally  (with  respect  to  M)  when  the  laser  is  moving  away  from  the  receiver.  The 
largest  value  of  TMAX  in  all  cases  occurs  when  4>0  =  180°  and  M  =2.  since  this 
combination  induces  the  greatest  stretching  of  the  received  signal. 

Figures  8-7(a)-(c)  represent  identical  cases  except  for  the  depth  of  the 
receiver.  They  show  that  T M  AX  is  sensitive  to  changes  in  h,  and  especially  when 
v  >  c. 

The  cases  used  to  generate  Figs.  8-7(e)-(g)  are  identical  except  for 
the  modulation  frequency  of  the  laser.  Like  TMIN ,  TMAX  decreases  as  the 
frequency  increases,  and  also  acts  almost  as  if  it  were  inversely  proportional  to  /. 

Figures  8-7(a),  (d),  and  (e)  were  generated  by  identical  situations  except 
for  the  source- to- receiver  range.  These  figures  show  that  TMAX  is  not  particu¬ 
larly  sensitive  to  changes  in  r0. 

From  these  results  we  conclude  that  the  time  property  TMAX  is  sensi 
five  to  changes  in  the  parameters  h,  /,  A/,  and  d>0.  Also,  like  T M IN,  the  manner 
in  which  TMAX  varies  with  respect  to  /,  M,  and  <^0  is  very  evident. 

8-1.6  AVERAGE  PERIOD,  TAVE 

Figure  8  8  shows  the  results  for  T AV E,  the  average  of  all  periods  in  the 
received  signal.  Like  the  two  previous  time  properties,  TAVE  is  very  sensitive 
to  changes  in  <p0  and  M.  Also,  TAVE  is  least  sensitive  to  changes  in  M  when 
0„  =  !>()°.  but  the  data  for  these  cases  indicate  that  TAVE  increases  slightly  as 
M  increases.  Like  TMIN  and  TMAX ,  TAVE  increases  monotonic  ally  (with 
respect  to  M)  when  the  laser  is  moving  away  from  the  receiver,  in  all  cases,  the 


largest  value  of  T AV E  occurs  when  4>q  =  180°  and  M  =2,  since  this  combination 
of  parameters  induces  the  greatest  elongation  of  the  received  signal. 

Figures  8-8(a) -(c)  were  generated  from  identical  cases  except  for  the 
depth  of  the  receiver.  These  figures  show  that  T  AV  E  is  sensitive  to  ft,  particularly 
when  A/  equals  1  or  1.5,  and  the  source  is  moving  towards  the  receiver. 

The  cases  used  to  generate  Figs.  8-8(e)-(g)  are  identical  except  for  the 
modulation  frequency  of  the  laser.  Like  TM I N  and  TM  AX ,  T AV E  decreases  as 
the  modulation  frequency  increases,  and  also  behaves  almost  as  if  it  were  inversely 
proportional  to  /. 

Figures  8-8(a),  (d),  and  (e)  represent  identical  situations  except,  for  the 
source- to- receiver  range.  These  figures  show  that  T AV E  is  not  very  sensitive  to 
changes  in  r0. 

From  these  results  we  conclude  that  T AV E  is  sensitive  to  changes  in 
the  parameters  ft,  /,  A/,  and  4>ot  and  that  its  behavior,  like  TMIN  and  TM  AX . 
with  respect  to  changes  in  /,  A/,  and  <f>0  is  very  evident. 

8-2  PRESSURE  RELATED  PROPERTIES 

The  pressure  related  properties  are  those  which  are  acquired  from  the 
pressure  peak  or  the  root  mean  square  pressure:  PPK.  P PM S.  and  PPK/PPMS. 

8-2.1  PRESSURE  PEAK,  PPK ,  AND  ROOT  MEAN  SQUARE  PRES¬ 
SURE,  PPMS 

Figure  8  f)  shows  t  he  results  for  PPK .  the  absolute  peak  pressure  of  the 
received  signal.  F  igures  8  !)(a)  (g)  show  t hat  PPK  is  very  sensitive  to  changes  in 
Od  and  M .  However,  notice  that  PPK  is  nol  sensitive  to  changes  in  On  when  the 


laser  is  stationary.  Also,  PPK  is  least  sensitive  to  changes  m  U  when  on 


')()'. 


since  at  this  point  the  sound  initially  appears  to  be  generated  by  a  stationary 
source.  Figures  8-9(a)-(g)  also  show  that  PPK  is  less  sensitive  to  changes  in 
M  when  the  laser  is  moving  away  from  the  receiver,  and  that  generally  PPK 
decreases  under  these  conditions  as  M  increases. 

Figures  8  9(a)  (c)  represent  identical  cases  except  for  the  depth  of  the 
receiver.  From  these  three  graphs  one  sees  that  PPK  generally  increases  as  li 
increases  when  M  equals  0.5  or  2,  and  decreases  with  increasing  h  when  M  equals 
1  or  1.5. 

The  cases  used  to  generate  Figs.  8  9(e)  (g)  are  identical  except  for 
the  modulation  frequency  of  the  laser.  These  figures  show  that  PPK  gener¬ 
ally  increases  as  /  increases.  This  result  can  be  explained  by  the  fact  that  the 
magnitude  of  the  time  derivative  of  the  laser  intensity  (which  drives  the  pressure) 
increases  as  /  increases. 

Figures  8  9(a),  (d).  and  (o)  were  generated  from  identical  situations 
except  for  the  source- to- receiver  range.  These  figures  show  that  PPK  generally 
decreases  as  r0  increases  due  to  spherical  spreading. 

From  these  results  we  conclude  that  the  absolute  peak  pressure  is  sensi¬ 
tive  to  changes  in  all  of  the  parameters. 

The  results  for  PPMS  are  shown  in  Fig.  8  10.  The  patterns  of  behavior 
as  parameters  are  varied  are  identical  for  PPMS  and  PPK .  Hence  we  conclude 
that  PPMS  is  sensitive  to  changes  in  all  the  parameters. 

8  2.2  PRESSURE  RATIO,  PPK/ PPMS 

I  lie  results  for  PPK  /PPMS  are  shown  in  Fig.  8  11.  Figures  S  11(a) 
(g)  show  that  PPh/PPM S  is  quite  sensitive  to  changes  in  o,,  and  M.  hike  all 
of  the  previous  piopert  ies.  PPK/PP.MS  is  not  sensitive  to  changes  in  the  angle 


do  when  the  laser  is  stationary.  Also,  notice  that  PPK/PRMS  is  sensitive  to 
changes  in  M  even  when  <^>0  =  90°. 

Figures  8--ll(a)-(c)  represent  identical  cases  except  for  the  depth  of  the 
receiver.  From  these  three  graphs  one  sees  that  PPK/PRMS  is  very  sensitive 
to  changes  in  h,  but  that  its  sensitivity  to  changes  in  h  decreases  as  h  increases. 

The  cases  used  to  generate  Figs.  8 -11(e) -(g)  are  identical  except  for  the 
modulation  frequency  of  the  laser.  These  figures  show  that  PPK/PRMS  is  quite 
sensitive  to  changes  in  /.  In  particular,  PPK/PRMS  seems  to  be  more  sensitive 
to  changes  in  /  when  the  laser  is  moving  towards  the  receiver. 

Figures  8  11(a),  (d),  and  (e)  were  generated  from  identical  situations 
except  for  the  source- to- receiver  range.  They  show  that  PPK/PRMS  is  very 
sensitive  to  changes  in  r0. 

From  these  results  we  conclude  that  the  ratio  PPK/PRMS  is  sensitive 
to  changes  in  all  the  parameters.  In  fact,  it  is  the  most  sensitive  measure  studied. 


CHAPTER  9 


SUMMARY  AND  CONCLUSIONS 


This  part  of  the  report  has  discussed  a  numerical  investigation  of  the 
generation  of  underwater  sound  by  a  moving  high  power  laser  source.  The  acoustic 
signal  generated  by  the  moving  laser  was  studied  in  detail.  Properties  of  the  signal 
were  identified  and  the  manner  in  which  these  properties  change  with  respect  to 
source  and  receiver  parameters  was  quantified.  This  information  identified  the 
aspects  of  the  sound  signal  that  are  least  and  most  susceptible  to  change. 

A  computer  program  called  MTS,  written  by  Berthelot,19  with  modi¬ 
fications  suitable  for  this  project,  was  used  to  numerically  predict  the  received 
acoustic  signal.  This  computer  program  is  based  on  the  time  domain  approach 
discussed  in  Chapter  1. 

The  parameters  of  the  source  which  were  varied  were  the  modulation 
frequency,  /,  of  the  laser  intensity  and  the  velocity,  t>,  or  Mach  number,  M  =  v/c, 
of  the  laser  beam  on  the  surface  of  the  water.  The  source- receiver  geometrical 
parameters  which  were  varied  were  the  initial  distance  between  the  source  and 
receiver.  r0,  the  initial  horizontal  angle,  <j)Q,  and  the  depth  of  the  receiver,  h. 
The  values  of  the  parameters  were  carefully  chosen  to  cover  the  whole  range  of 
physically  realizable  situations.  All  possible  combinations  of  the  parameters  were 
used,  resulting  in  a  total  of  67.r)  cases. 
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Ihe  properties  of  the  received  acoustic  signal  which  were  examined  may 
he  placed  into  two  categories:  time  related  and  pressure  related  properties.  The 
time  related  properties  of  the  acoustic  signal  which  were  numerically  investigated 
were  the  duration  of  the  received  signal  TD.  the  time  at  which  the  absolute  peak 
pressure  occurs  7  PA  .  the  maximum  period  in  the  signal  TM  AX.  the  minimum 
period  /  \//.\  ,  the  average  of  all  periods  TA'VE,  and  the  time  inversion  prop¬ 
erty  TPh  /I  D.  I  he  pressure  related  properties  studied  were  the  absolute  peak 
pressure  /’PA',  the  root  mean  square  pressure  of  the  acoustic  signal  PP.\/ S.  and 
the  pressure  ratio  PPA/P RMS. 

I  he  results  showed  that  TD  is  very  sensitive  to  changes  in  the  initial 
horizontal  angle,  <?0,  and  the  Mach  number,  M.  of  the  laser  beam  on  the  surface 
of  the  water.  However,  77)  is  not  sensitive  to  changes  in  the  angle  0O  when  the 
source  is  stationary.  Also.  TD  is  least  sensitive  to  changes  in  M  when  C*o  =  !)0o.  In 
the  case  when  0o  =  9Q°.  the  Mach  number  as  seen  by  the  receiver.  RM .  is  initially 
equal  to  zero,  and  remains  approximately  equal  to  zero  for  the  duration  of  the 
laser  pulse.  To  the  receiver,  then,  the  sound  essentially  appears  to  be  generated 
by  a  stationary  source.  Hence  the  duration  of  the  received  pulse  is  relatively 
constant . 

It  was  also  found  that  TD  increases  monotonically  (with  respect  to  M) 
as  the  laser  moves  away  from  the  receiver.  The  maximum  duration  of  the  received 
signal  occurs  when  the  laser  is  moving  directly  away  from  the  receiver.  The  results 
also  showed  that  TD  is  most  sensitive  to  changes  in  M  when  the  laser  is  moving 
directly  towards  the  receiver,  or  directly  away  from  the  receiver.  However,  the 
results  showed  that  I  D  is  very  insensitive  to  changes  in  the  parameters  r0,  f. 
and  h. 


from  these  results  it  was  concluded  that  TD  is  dependent  on  the  apparent 
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speed  of  the  source.  In  terms  of  geometrical  parameters,  this  means  that  the 
duration  of  the  received  signal  is  highly  dependent  on  0O  and  M . 

Like  the  time  property  ID,  the  results  showed  that  TPK  is  quite  sensi¬ 
tive  to  changes  in  <p0  and  A/.  For  a  stationary  source,  TPK .  like  TD,  is  not 
a  function  of  g>0.  When  the  laser  is  moving  directly  towards  the  receiver,  the 
results  showed  that  TPK  decreases  as  RM  goes  from  0  to  1,  and  increases  as 
RM  goes  from  1  to  2.  This  is  a  result  of  the  signal  being  compressed  in  time, 
as  RM  —  1  is  approached  from  below  and  above.  It  was  found  that  TPK  is  not 
symmetric  about  RM  =  1.  A  second  effect,  resulting  from  inversion  of  the  signal 


when  RM  >  1.  skew’s  the  curve. 


The  results  also  showed  that  TPK  is  sensitive  to  changes  in  /<,  especially 
when  the  velocity  of  the  laser  on  the  surface  of  the  water  is  greater  than  the  speed 
of  sound  in  water.  It  was  also  found  that  TPK  is  sensitive  to  changes  in  f.  even 
when  M  =0.  Except  when  A/  =  0,  T PK  is  sensitive  to  changes  in  r0,  and  more  so 
when  v>c.  From  these  results  it  was  concluded  that  TPK  is  sensitive  to  changes 
in  all  the  parameters. 

The  results  showed  that  the  time  inversion  property,  T PI\ /T D,  is  gener¬ 
ally  not  sensitive  to  changes  in  4>o  when  v  <  c,  but  is  quite  sensitive  to  changes  in 
Ox  when  v  >c,  and  the  laser  is  moving  towards  the  receiver.  Also,  TPK/TD  is 
generally  not  sensitive  to  changes  in  M  when  the  laser  is  moving  away  from  the 
receiver,  but  is  senutive  to  changes  in  M  when  the  laser  is  moving  towards  the 
receiver.  I  he  results  also  showed  that  TPK/TD  is  fairly  insensitive  to  changes 
in  li,  f,  and  r„.  From  these  results  it  was  concluded  that  the  received  signal  is 


heavily  lime  inverted  when  tin’  laser  velocity  is  greater  than  the  sound  speed,  and 


is  moving  towards  the  receiver.  In  addition,  the-  level  of  time  inversion  is  highly 
dependent  on  o, i  and  M. 
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The  property  TMIN  was  found  to  be  quite  sensitive  to  changes  in  <f>0  and 
M.  However,  T MIN  is  not  sensitive  to  changes  in  <f>0  when  the  laser  is  stationary. 
Also,  no  change  in  TMIN  was  detected  as  M  was  varied  when  <^o  =  90°.  Again, 


> 
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for  this  case  the  signal  appears  to  be  generated  by  an  almost  stationary  source. 
The  results  also  showed  that  TMIN  is  particularly  sensitive  to  changes  in  the 
velocity  of  the  laser  when  it  is  moving  directly  towards  and  away  from  the  receiver. 

It  was  also  found  that  TMIN  increases  monotonically  (with  respect 
to  M)  when  the  laser  is  moving  away  from  the  receiver.  The  largest  value  of 
TMIN  occurs  when  the  laser  is  moving  directly  away  from  the  receiver,  since 
this  condition  produces  the  greatest  stretching  of  the  received  signal. 

The  results  showed  that  TMIN  is  sensitive  to  changes  in  h,  particularly 
when  the  source  is  moving  towards  the  receiver,  and  M  equals  1  or  1.5.  As 
expected,  TMIN  was  found  to  act  almost  as  if  it  were  inversely  proportional  to 
/,  regardless  of  the  Mach  number.  The  results  also  showed  that  TMIN  is  not 
particularly  sensitive  to  changes  in  r0. 

The  patterns  of  behavior  as  parameters  were  varied  for  T M AX  and 
TAVE  were  found  to  be  practically  identical  to  those  of  TMIN.  The  only 
difference  in  their  behavior  was  found  to  occur  when  <f>o  =  90°.  When  d>o  =  90°, 
both  T M AX  and  TAVE  were  found  to  increase  slightly  as  M  was  increased. 

From  these  results,  it  was  concluded  that  TM AX  and  TAVE  are  sensi¬ 
tive  to  changes  in  the  parameters  /«,  /,  A/,  and  <t> o-  Furthermore,  the  manner  in 
which  I'M  AX  and  TAVE  vary  with  respect  to  /,  A/,  and  4>0  is  very  evident. 

The  numerical  results  showed  that  P P I\  is  very  sensitive  to  changes  in 
C>o-  However,  PPK  was  found  to  be  insensitive  to  changes  in  <j f>0  when  the  laser 
is  stationary,  since  at  this  point  the  sound  initially  appears  to  be  generated  by  a 
stationary  source.  Also,  PPK  is  less  sensitive  to  changes  in  M  when  the  laser  is 
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moving  away  from  the  receiver,  and  generally  decreases  under  these  conditions. 

The  results  also  showed  that  PPK  generally  increases  as  h  increases  and 
as  /  increases.  The  latter  can  be  explained  by  the  fact  that  the  magnitude  of  the 
time  derivative  of  the  laser  intensity  (which  drives  the  pressure)  increases  as  / 
increases.  Furthermore,  the  results  showed  that  PPI\  generally  decreases  as  7-0 
increases.  This  decrease  is  due  to  spherical  spreading.  From  these  results  it  was 
concluded  that  PPI\  is  sensitive  to  changes  in  all  of  the  parameters. 

The  results  for  PPM S  showed  that  its  patterns  of  behavior  as  parame¬ 


ters  are  varied  are  identical  to  those  of  PPK.  It  was  concluded  that  PRMS  is 


sensitive  to  changes  in  all  the  parameters. 

Finally,  the  results  showed  that  PPK / PRMS  is  very  sensitive  to  changes 
in  do  and  M.  Like  all  of  the  previous  properties,  PPK/ PRMS  is  not  sensitive 
to  changes  in  the  angle  do  when  the  laser  is  stationary.  Also,  PPK/ PRMS  was 
found  to  be  sensitive  to  changes  in  M  even  when  do  =  90°.  The  results  also  showed 
that  PPK/ PRMS  is  very  sensitive  to  changes  in  h  and  r0.  Furthermore,  it  was 
found  that  PPK/ PRMS  is  quite  sensitive  to  changes  in  /,  especially  when  the 
laser  is  moving  towards  the  receiver.  From  these  results  it  was  concluded  that  the 
ratio  PPK/ PRMS  is  sensitive  to  changes  in  all  the  parameters,  and  that  this 
measure  is  more  sensitive  than  any  other  studied. 

Taken  as  a  whole,  this  numerical  study  has  shown  that  <ji>o  and  M  are  the 
most  critical  parameters,  affecting  virtually  all  of  the  received  signal  properties. 
In  addition,  f  is  a  critical  parameter  for  all  of  the  pressure  related  properties,  and 
the  time  related  properties  derived  from  the  zero  crossing  of  the  received  signal. 
Hence,  if  thermoacoustic  generation  of  sound  is  to  be  used  for  communication, 
extreme  care  must  be  taken  that  M ,  do,  and  f  are  well  known. 
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